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EXECUTIVE  SUMMARY 

A.  INTRODUCTION 

This  paper  was  prepared  for  the  Federal  Aviation  Adminis¬ 
tration,  High  Altitude  Pollution  Program  as  a  part  of  the  con¬ 
tracted  task.  Analysis  of  Aircraft  Effects:  Natural  Atmospheric 
Perturbation. 

Present  and  projected  transport  aircraft  operations  will 
inject  a  substantial  amount  of  oxides  of  nitrogen  (N0X  *  NO 
+  NO2)  into  the  atmosphere  at  cruise  altitudes.  The  signifi¬ 
cance  of  these  injections  on  atmospheric  ozone  and  overall 
atmospheric  chemistry  can  only  be  determined  by  comparing  the 
aircraft  source  with  all  other  sources,  natural  and  anthropo¬ 
genic.  A  comparison  of  the  N0X  source  strengths  shown  in  Table 
S— 1  indicates  that  lightning  is  a  significant  natural  source  of 
tropospheric  N0X.  This  paper  quantifies  the  injection  rate  of 
N0X  due  to  lightning  and  its  variation  with  altitude,  geographic 
location,  and  season,  pointing  out  significant  uncertainties 
in  the  estimates. 

A  global  lightning  rate  of  300  flashes/s  has  been  adopted 
here,  consisting  of  approximately  20  percent  cloud-to-ground 
flashes  and  80  percent  intracloud  discharges.  A  production 
rate  of  1026  NOx/flash  for  ground  flashes  (Noxon  1976;  1978), 
and  1025  NOx/flash  for  cloud  discharges  yields  an  annual  global 
production  rate  of  5.7  Tg  N/yr  ( 1  Tg  =  1012  grams)  from  lightning 
with  an  estimated  uncertainty  range  of  2-20  Tg  N/yr. 

The  discussion  begins  with  a  description  of  four  forms  of 
atmospheric  electricity  related  to  lightning,  listed  here  in 
the  order  of  their  energy  dissipation  per  unit  length;  corona 
discharges,  intracloud  discharges,  cloud-to-ground  strokes  and, 
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superbolt  flashes.  The  net  N0X  production  by  electrical  dis¬ 
charges  is  proportional  to  the  discharge  energy.  Corona  dis¬ 
charges  ,  the  low-energy  discharges  which  can  occur  almost  con¬ 
tinuously  during  a  thunderstorm,  actually  contribute  only  mini¬ 
mally  to  the  N0X  production.  At  the  other  extreme,  superbolt 
flashes  are  extremely  energetic  but  rare,  so  that  their  cumula¬ 
tive  N0X  production  rate  is  also  minimal.  The  rest  of  the 
report  concerns  itself  with  the  N0X  production  from  intracloud 
discharges  and  cloud-to-ground  flashes. 

TABLE  S-l.  SOURCE  STRENGTHS  OF  TROPOSPHERIC  NOx* 

Tg  N/yr 

Natural  Sources: 


Lightning  5.7 

Forest  Fires  1.7 

Transport  from  Stratosphere  0.5 

Cosmic  rays  0.06 

Exhalation  from  soils  10 

Anthropogenic  Sources  (1975  Estimates) : 

Fossil  fuel  combustion  20 

Biomass  burning**  3.3 

Aircraft  0.15 

NH3  decomposition  (also  natural)  <8 


•Compiled  from  Bauer  (1982) 

••includes  deforestation  due  to  population  increase, 
shifting  agriculture,  and  the  burning  of  firewood, 
bush,  and  savannah. 

B.  NOx  PRODUCTION  BY  LIGHTNING 

The  amount  of  NO)t  produced  by  lightning  has  been  measured 
in-situ,  simulated  in  the  laboratory  with  spark  discharges,  and 
numerically  modeled  by  various  methods.  Most  of  the  final  esti 
mates  are  between  (0.6-4)  x  10^  NOx  molecules  produced  per 
lightning  flash. 

This  study  adopts  the  recommendation  from  three  in-situ 
measurements  by  Noxon  (1976,  1978).  Under  the  dark  sky  condi¬ 
tions  during  a  thunderstorm  the  NO  produced  by  lightning  is 
largely  converted  to  NO2;  in  this  study  NO2  is  assumed  to  be 
equivalent  to  NOx.  Noxon  recorded  the  absorption  spectrum  of 
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the  overhead  column  of  NO2  before,  during,  and  after  a  thunder¬ 
storm.  Noxon's  (1976)  first  measurement  of  10^6  N02/flash  was 
supported  by  two  subsequent  experiments  (Noxon  1978),  and  is  in 
reasonable  agreement  with  most  of  the  numerical  model  simula¬ 
tions  of  lightning  in  air.  In-situ  measurements  such  as  Noxon’s 
have  several  advantages  for  an  inventory.  First,  the  lightning 
process  is  genuine;  therefore,  assumptions  about  lightning 
parameters  such  as  energy  deposited  per  stroke  and  additional 
NOx  produced  by  subsequent  strokes  in  a  flash  need  not  be 
made.  Second,  the  problem  of  rainout  immediately  after  the 
flash  need  not  be  considered,  because  the  in-situ  measurements 
give  a  net  NOx  production. 

C.  GLOBAL  DISTRIBUTION  OF  N0X  FROM  LIGHTNING 

In  this  work,  global  distributions  of  lightning  obtained 
from  the  DMSP  and  ISS-b  satellites*  have  been  emphasized,  al¬ 
though  in  general,  an  excellent  qualitative  agreement  exists 
with  the  historical  distribution  of  thunderstorms.  It  shoulc. 
be  noted  that  the  ISS-b  satellite  measurement  by  Kotaki  et  al. 
(1981)  was  done  by  sensing  the  radio-frequency  emissions  from 
lightning  discharges,  which  is  an  experimental  method  sensitive 
to  both  cloud-to-ground  flashes  and  intracloud  discharges, 
whereas  the  DMSP  satellite  measurement  by  Turman  and  Edgar 
(1980)  was  done  with  an  optical  sensor,  a  method  which  may 
discriminate  against  the  less  energetic  intrecloud  discharges. 

The  globally  integrated  rate  of  lightning  obtained  from 
these  two  measurements  differ  by  a  factor  of  three,  i.e.,  DMSP 
=80+40  flashes/s  and  ISS-b  =  280  flashes/s.  The  discrepancy 
between  these  two  determinations  may  be  due  to  the  normalization 
required  to  convert  the  DMSp  satellite  data  to  an  absolute 
lightning  flash  rate.  Specifics  about  the  normalization  are 
discussed  in  the  main  text. 

*DMSP  =  Defense  Meteorological  Satellite  Program,  USA 
ISS-b  =  Ionospheric  Sounding  Satellite,  Japan. 
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D.  INTRACLOUD  NOx  PRODUCTION  EFFICIENCY 

The  efficiency  of  NOx  formed  by  lightning,  as  well  as  its 
altitude  of  injection  in  the  atmosphere,  differs  for  cloud  dis¬ 
charges  and  ground  strokes.  In  this  work,  N0X  production  by 
intracloud  discharges  was  assumed  to  be  one-tenth  as  efficient 
when  compared  to  a  ground  flash,  based  on  experimental  evidence 
that  the  amount  of  N0X  produced  is  proportional  to  the  discharge 
energy.  This  assumption  is  a  major  source  of  uncertainty,  since 
no  direct  experimental  verification  of  the  energy  deposited  by 
intracloud  discharges  exists. 

E.  INJECTION  HEIGHT  AND  RESIDENCE  TIME 

The  altitude  of  N0X  production  is  not  the  altitude  of  N0X 
injection  because  of  the  violent  updrafts  and  downdrafts  asso¬ 
ciated  with  the  center  of  a  thunderstorm.  The  stabilization 
height  was  taken  to  be  directly  below  the  local  tropopause, 
that  is,  12  km  for  midlatitudes  and  15  km  for  the  tropics.  A 
density-weighted  vertical  distribution  was  proposed  for  the 
N0X  injected  by  both  intracloud  and  cloud-to-ground  light¬ 
ning.  Experimental  evidence  in  which  the  luminosity  of  the 
lightning  channel  decreased  with  height  and  the  tendency  to 
produce  a  density-weighted  vertical  distribution  with  suffi¬ 
cient  atmospheric  mixing  support  this  assumption.  The  N0X 
from  intracloud  discharges  was  assumed  to  be  distributed  in  a 
5-km  band  below  the  local  tropopause  and  the  N0X  from  ground 
flashes  was  distributed  uniformly  between  that  level  and  the 
ground. 

The  effect  o 1  chemical  injections  on  atmospheric  photo¬ 
chemistry  depends  on  the  residence  time  of  the  chemical  species, 
which,  in  turn,  depends  on  the  altitude  and  latitude  of  the 
injection.  N0X  injected  at  high  tropospheric  altitudes,  espe¬ 
cially  near  the  tropopause,  will  remain  in  the  atmosphere  Signif¬ 
icantly  longer  than  an  equal  amount  from  a  source  near  the  ground. 
Similarly,  the  residence-  time  of  N0X  in  regions  where  the  net  flow 


in  atmospheric  circulation  is  upwards,  as  in  the  tropics,  will 
be  longer  than  in  regions  of  net  downward  flow. 

F.  INTERNAL  RAINOUT 

Various  mechanisms  for  N0X  loss  at  ambient  temperatures 
were  examined  and  were  found  to  be  inefficient  on  a  time  scale 
comparable  to  the  lifetime  of  a  thunderstorm.  The  NO2  produc¬ 
tion  estimate  from  Noxon's  {1976;  1978)  in-situ  measurements 
can  be  considered  to  be  the  residual  NO2  which  remains  after 
any  initial  rainout  of  HNO3  or  HNO2  formed  in  the  hot  liqhtning 
channel.  Any  losses  during  the  thunderstorm  that  occur  after 
Noxon's  measurements  represent  true  losses  to  the  inventory  but, 
as  mentioned  previously,  these  losses  were  found  to  be  small. 

G.  DATA  FOR  2-D  MODEL 

The  recommended  treatment  for  the  2-D  distribution  is  sum¬ 
marized  here.  An  annual  average  lightning  rate  of  300  flashes/s, 
including  both  intracloud  discharges  and  cloud- to-ground  flashes, 
was  chosen.  This  rate  was  partitioned  among  the  latitudes  and 
months  of  the  year  according  to  the  global  lightning  distribu¬ 
tions  observed  by  Turman  and  Edgar  (1980)  from  the  DMSP  satellite. 

The  monthly  values  can  be  appropriately  averaged  for  sea¬ 
sonal  or  annual  calculations .  The  values  at  each  latitude  were 
in  turn  partitioned  between  ground  flashes  and  intracloud  dis¬ 
charges  based  on  the  recommendations  of  Prentice  and  Mackerras 
(1977).  The  upper  bound  of  the  intracloud  injection  region  is 
defined  by  the  local  tropopause.  For  use  in  a  numerical  model, 
the  NOx  production  rate  in  each  injection  region  must  be 
distributed  among  the  number  of  altitude  bins  available  in  the 
model.  A  density-weighted  distribution  is  recommended.  Table 
S-2  specifies  the  annual  NOx  injection  rate  in  Tg  N/yr  as  a 
function  of  latitude  and  altitude. 

Several  uncertainties  can  be  identified  in  the  estimate 
of  the  lightning  contribution  to  the  global  NO.,  budget  in  the 


TABLE  S-2 .  ANNUAL  N0X  INJECTION  RATE  FROM 
LIGHTNING:  2-D  DISTRIBUTION 


1()12  gN/yr-tan* 


Atmospheric 
Nunber 
Density** 
x  10«  km-3 


0.045 

0.033 

0.030 

0.057 

0.039 

0.035 

0.060 

0.045 

0.041 

0.072 

0.053 

0.049 

0.081 

0.060 

0.055 

0.030 

0.024 

0.024 

0.034 

0.027 

0.027 

0.038 

0.031 

0.030 

0.043 

0.034 

0.034 

0.048 

0.038 

0.038 

0.053 

0.043 

0.042 

0.059 

0.048 

0.046 

0.065 

0.053 

0.051 

0.072 

0.058 

0.057 

0.079 

0.064 

0.063 

0.032 

0.012 

0.001 

0.037 

0.014 

0.002 

0.042 

0.016 

0.002 

0.047 

0.018 

0.002 

0.053 

0.020 

0.003 

0.060 

0.034 

0.005 

0.067 

0.038 

0.006 

0.075 

0.042 

0.006 

0.083 

0.046 

0.007 

0.092 

0.051 

0.007 

0.101 

0.057 

0.008 

0.112 

0.063 

0.009 

LATITUDE 


The  aonal  areas  within  XO-degree  latitude  intervals  are  given  in  Table  10  to  facilitate 
conversion  to  injection  units  such  as  molecules/aiPs. 


U.S.  Standard  Atmosphere,  1976. 


troposphere: 

(1)  The  N0X  production  efficiency  for  intracloud 
discharges 

(2)  The  N0X  production  efficiency  in  the  tropics 

(3)  The  injection  height  of  N0X  produced  during  a  light¬ 
ning  storm  both  at  midlatitudes  and  in  the  tropics 

(4)  The  normalization  of  the  global  lightning  rate  meas¬ 
ured  by  satellites,  especially  in  the  tropics. 

(5)  The  ratio  of  cloud  discharges  to  ground  flashes,  also 
with  particular  emphasis  in  the  tropics  where  current 
data  is  scarce. 

The  accumulated  uncertainty  for  the  measured  quantities 
has  been  estimated  to  be  within  one  order  of  magnitude.  It  is 
not  possible  to  estimate  the  uncertainty  for  the  unmeasured 
quantities  such  as  the  N0X  production  efficiency  by  intracloud 
discharges  and  the  injection  height. 

H.  RECOMMENDATIONS 

Several  experiments  that  could  eliminate  most  of  the  un¬ 
certainty  in  the  lightning-generated  N0X  estimate  are  recom¬ 
mended.  In  most  cases,  these  can  be  carried  out  with  currently 
developed  equipment  and  some  minor  modifications. 

(1)  The  F-106B  aircraft  used  in  the  NASA-Langley  Storm 
Hazards  Project,  which  is  currently  equipped  to  fly 
within  thunderstorms,  should  be  fitted  to  measure  N0X 
at  various  altitudes  within  the  storm  cloud,  at  the 
perimeter  of  the  cloud,  and  outside  the  cloud.  This 
would  answer  the  questions  concerning  N0X  production 
efficiency  as  well  as  those  concerning  the  height  of 
injection  in  the  region  of  net  updraft  in  the  center 
of  the  thundercloud  and  in  the  region  of  net  downdraft 
at  the  perimeter.  Continued  measurements  as  the  storm' 
passes  would  indicate  the  rate  of  N0X  dissipation  and 
chemical  removal. 
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(2)  The  tropical  regions  should  be  included  if  any  further 
experiments  of  the  N0X  production  efficiency  by  light¬ 
ning,  such  as  Noxon's  column  measurements,  are  done. 

(3)  The  calibration  of  lightning  flash  rate  measured  by 
satellites  needs  more  serious  consideration.  The 
several  global  maps  generated  by  satellite  have  simi¬ 
lar  distributions,  but  a  consensus  for  the  globally 
averaged  value  does  not  exist.  Any  ground-based 
calibration  data  must  include  tropical  as  well  as  mid¬ 
latitude  measurements.  Attention  must  also  be  given 

to  the  ratio  of  intracloud  discharges  to  ground  flashes, 
Nc/Ng  for  the  calibration  points  chosen. 

(4)  The  uncertainty  in  the  ratio  of  the  number  of  intracloud 
discharges/ground  flashes,  Nc/Ng,  for  tropical  regions 
needs  to  be  reduced  by  additional  measurements.  It  is 
emphasized  that  most  lightning  occurs  in  the  tropics 
but  current  measurements  are  concentrated  in  the  mid¬ 
latitudes. 
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I.  INTRODUCTION 

To  understand  the  significance  of  anthropogenic  NOx  emis¬ 
sions,  the  magnitude  and  altitude  of  the  emissions  must  be 
compared  to  natural  sources  within  the  atmosphere.  The  con¬ 
tribution  of  lightning  has,  to  date,  been  a  major  uncertainty 
in  such  a  global  N0X  source  inventory.  Current  estimates  of 
the  N0X  production  rate  by  lightning  vary  by  a  factor  of  50 
from  1.8  to  90  Tg  N/yr.  This  range,  however,  is  not  indicative 
of  the  uncertainty  of  global  NOx  production  by  lightning  since 
not  all  the  estimates  were  made  with  the  same  initial  conditions 
and  assumptions  about  the  lightning  process.  If  only  cloud-to- 
ground  flashes  are  counted,  the  estimates  are  between  1.8  and 
22  Tg/yr  and  if  intracloud  lightning  is  included,  the  range  be¬ 
comes  9  to  90  Tg/yr.  A  comparison  with  anthropogenic  combustion 
sources,  estimated  at  about  20  Tg  N/yr  (Crutzen  1979),  shows 
that  these  estimates  span  the  entire  range  from  a  minor  to  a 
significant  source  of  NOx.  This  paper  examines  the  separate 
components  necessary  to  make  an  estimate  of  the  global  production 
rate  of  NOx  by  lightning.  Previous  work  in  which  global  NOx 
production  rates  were  estimated  aire  examined,  with  particular 
emphasis  on  the  experimental  and  theoretical  basis  of  previous 
estimates.  Recommendations  are  made  based  on  the  least  ambiguous 
data  source.  The  remainder  of  this  chapter  will  discuss  some 
basic  information  about  the  lightning  process  and  other  forms 
of  atmospheric  electricity  that  may  be  sources  of  NOx  in  the 
troposphere. 

It  is  clear  to  anyone  who  has  seen  the  damage  to  an  object 
struck  by  lightning  that  massive  energies  are  released.  During 
a  storm,  large  potential  energies  are  created  when  regions  of 
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positive  and  negative  charge  accumulate  within  a  cloud.  The 
exact  mechanism  of  the  charge  generation,  separation,  and 
accumulation  is  unknown,  but  most  theories  consider  the  process 
to  be  due  either  to  the  frictional  effects  of  falling  water  and 
ice  particles  within  a  cloud  or  to  the  transport  and  accumula¬ 
tion  of  external  charge  by  convective  circulation  within  the 
cloud.  When  sufficient  charge  has  accumulated,  on  the  order 
of  hundreds  of  coulombs  within  a  thundercloud  a  few  kilometers 
wide,  the  high  electric  field  and  space-charge  effects  can 
result  in  a  point,  or  corona,  discharge.  An  arc  discharge  is 
created  if  a  sufficient  potential  difference  exists  between 
the  cloud  and  the  far  edge  of  the  corona  discharge.  A  portion 
of  the  cloud's  charge  can  then  be  transferred  to  the  leading 
point  of  the  newly  created  ionized  path.  The  process  can  be 
repeated  in  a  stepwise  manner  so  that  an  atmospheric  breakdown 
path  many  kilometers  long  can  result,  even  though  a  direct  one- 
step  breakdown  would  not  be  possible.  The  partially  ionized 
breakdown  channel  is  called  the  stepped  leader. 

A  potential  for  massive  charge  neutralization  exists  if  the 
stepped  leader  reaches  the  ground.  The  neutralization  that 
begins  at  the  ground  and  then  progresses  up  the  stepped  leader 
path  is  referred  to  as  the  return  stroke.  About  90  percent  of 
the  lightning  stroke  energy  is  dissipated  in  the  return  stroke, 
but  only  a  fraction  of  the  cloud's  total  charge  is  neutralized. 
Subsequent  strokes  along  the  same  lightning  channel  can  neutral¬ 
ize  the  residual  charge. 

A  typical  number  of  strokes  per  lightning  flash  is  three 
to  four;  however,  the  number  depends  on  factors  such  as  light¬ 
ning  type,  storm  type,  and  specific  geographic  location.  A 
clear  dependence  on  latitude  is  not  evident.  The  stroke-per- 
flash  data  summarized  in  Table  1  shows  that  multiple  strokes 
are  prevalent  in  New  Mexico  and  single  flashes  are  more  common 
on  Mount  San  Salvatore,  Switzerland.  Schonland  (1956)  reported 
that  50  percent  of  the  flashes  out  of  a  total  of  1800  in  South 
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Africa  had  four  or  more  strokes  per  flash.  He  also  noted  that 
multiple  strokes  occur  more  readily  in  large  frontal  storms 
than  in  small  convective  storms. 


TABLE  1.  NUMBER  OF  STROKES  WITHIN  A  FLASH 


Number  of 
Strokes 
per  Flash 

Percentage  of  Total  Flashes 

South  Africa3 

New  Mexico*3 

Switzerland0 

1 

13 

12.2 

76.5 

2 

19 

26.5 

7.7 

3 

18 

7.3 

4.8 

4 

20 

4.9 

1.9 

5 

12 

- 

2.3 

6 

14.6 

1.4 

7 

18 

2.5 

1.0 

8 

4.9 

0.7 

>8 

27. 1 

3.8 

Total  Number 

of  flashes 

530 

41 

1026 

Mean  value: 

strokes/flash 

1220/530*2.3 

265/41*6.5 

1950/1026=1.9 

a  Malan  (1956),  cited  in  Berger  (1977). 

b  Kitagawa  et  al.  (1962),  cited 

in  Berger  (1977) . 

c  Berger  (1972),  cited  in  Berger  (1977). 

An  intracloud  discharge  can  occur  if  a  leader  within  a 
cloud,  called  a  streamer,  reaches  a  point  with  a  significantly 
lower  potential,  i.e.  ,  a  differently  charged  portion  within  the 
same  cloud,  or  to  another  cloud.  The  current  of  an  intracloud 
discharge,  about  100  amperes,  is  significantly  less  than  the 
104  to  105  amperes  characteristic  of  the  cloud-to-ground 
return  stroke.  It  should  be  noted  that  cloud-to-ground  light¬ 
ning  generally  has  an  intracloud  component  but  not  all  intracloud 
lightning  eventually  hits  the  ground. 


Most  corona  discharges  do  not  develop  into  cloud-to-ground 
lightning  or  intracloud  and  intercloud  discharges.  Typically, 
the  energy  dissipated  by  a  corona  discharge  is  low,  with  currents 
on  the  order  of  microamperes.  Berger  (1977)  cites  an  example 
of  an  experiment  on  San  Salvatore  where  corona  discharges  trans¬ 
ported  a  total  of  72  coulombs  to  the  experimental  tower  within 
20  hours.  This  is  comparable  to  the  charge  transfer  within  a 
single  lightning  stroke.  It  must  be  noted  that  the  corona  dis¬ 
charge  process  occurs  continuously  over  a  large  surface  area, 
therefore  the  time  and  space  integrated  charge  transfer  can  be 
quite  large.  For  example,  the  charge  transferred  from  the 
clouds  to  the  ground  by  corona  discharges  has  been  estimated 
from  an  integration  of  current  records  and  was  found  to  be  on 
the  average  20  times  greater  than  the  charge  transferred  by 
cloud-to-ground  lightning  strokes  (Latham  and  Stromberg  1977). 
Even  weak  electrical  discharges  can  cause  molecular  dissocia¬ 
tion  and  recombination,  as  is  evident  from  experimental  dis¬ 
charges  in  primitive  atmospheres  in  which  amino  acids  are 
synthesized  over  long  periods  of  time. 

Routine  satellite  coverage  of  the  earth  has  revealed  the 
existence  of  a  rare  type  of  lightning  called  a  superbolt  flash 
in  which  the  energy  dissipated  is  over  100  times  more  powerful 
than  typical  lightning  (Turman  1977).  For  over  5  years  the 
sensors  aboard  the  Vela  V  satellites  recorded  an  average  of 
5  superbolt  flashes/year  with  powers  greater  than  about 
3  x  10^2  watts  (Turman  1979).  This  low  lightning  rate  is  in¬ 
sufficient  to  affect  the  global  NOx  production,  however  it 
may  cause  local  enhancement  near  Japan  where  most  superbolts 
occur.  The  existence  of  superbolts  bring  out  a  crucial  but 
unanswered  question;  are  the  superbolts  a  separate  physical 
phenomenon,  or  are  they  the  high-energy  part  of  an  energy  dis¬ 
tribution  of  lightning?  If  the  former  is  true,  the  superbolts 
have  a  negligible  effect  on  the  annual  N0X  production  by  light¬ 
ning.  If  the  latter  is  true,  then  the  wide  energy  distribution 
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makes  it  more  difficult  to  reach  conclusions  about  a  "typical” 
lightning  flash.  In  this  case  one  must  be  careful  to  work  with 
an  adequate  statistical  sample. 

NOx  (NOx  =  NO  +  NO2)  production  can  occur  for  most  types 
of  atmospheric  discharges;  the  amount  of  N0X  production  depends 
on  the  discharge  energy  and  the  duration  of  the  discharge.  To 
date,  the  studies  of  N0X  production  by  lightning  have  concen¬ 
trated  on  the  cloud-to-ground  lightning  stroke.  On  energetic 
grounds  alone  this  approach  may  be  justified,  since  the  powerful 
return  stroke  of  cloud- to-ground  lightning  releases  the  most 
energy  and  produces  the  most  N0X  per  unit  distance.  Experimental 
discharge  studies,  discussed  in  Chapter  II,  indicate  that  the  N0X 
yield  per  joule  is  approximately  constant  for  low  and  high  energy 
discharges  between  10  and  10^  J/m,  which  implies  that  the  NOx 
produced  is  directly  proportional  to  the  discharge  energy.  Al¬ 
though  corona  discharges  occur  continuously  during  a  lightning 
storm  the  discharge  energy  is  quite  small.  Measurements  by 
Noxon  (1978)  indicate  that  corona  discharges  are  inefficient 
at  NOx  production.  The  measured  NO2  concentrations  during 
corona  activity  were  not  too  different  from  typical  ambient 
values . 

In  order  to  appreciate  the  magnitude  of  NOx  production 
contributed  by  each  type  of  discharge,  it  is  necessary  to  con¬ 
sider  the  cumulative  effects,  discharge  times,  distances,  and 
frequencies.  Although  intracloud  discharges  contain  no  large 
energy  sink  like  the  return  stroke  they  can  travel  over  longer 
distances.  Intercloud  strokes  in  excess  of  50  kilometers 
have  been  recorded.  It  is  also  necessary  to  consider  the 
frequency  of  discharge  events.  Intracloud  discharges  occur 
up  to  six  times  as  often  as  cloud-to-ground  strokes  and,  as 
mentioned  previously,  corona  discharges  occur  almost  continuously 
during  a  thunderstorm.  This  will  be  considered  in  more  detail 
in  Chapter  3  when  lightning  counting  is  discussed.  Several 
books  have  been  written  recently  about  the  basic  properties  of 
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II.  NOx  PRODUCTION  BY  LIGHTNING 


Estimates  of  the  amount  of  N0X  produced  per  unit  energy 
by  a  lightning  stroke  have  been  obtained  by  numerical  simula¬ 
tion  with  various  model  approaches,  laboratory  simulations 
with  spark  discharges  in  air,  and  in  one  case  with  an  in-situ 
measurement  of  NO2  in  a  lightning  storm.  Table  2  shows  that 
the  range  of  these  estimates  span  an  order  of  magnitude.  This 
section  will  briefly  discuss  the  techniques  used  to  estimate 
NOx  production  specifically  addressing  the  variability  evident 
in  Table  2. 

A.  NUMERICAL  MODEL  CALCULATIONS 
1 .  A  Description  of  Lightning  Models 

NO  production  by  lightning  is  modelled  by  two  different 
approaches:  the  shock  heating  model  proposed  by  Tuck  (1976) 

and  extended  by  Chameides  et  al.  (1977)  and  Chameides  (1979a, 
1979b)  and  the  hot  channel  model  proposed  by  Hill  et  al.  (1980) 
Both  models  calculate  NOx  production  by  the  number  of  NO  mole¬ 
cules  produced  in  air  at  high  temperatures  which  are  then 
frozen  in  as  the  air  expands  and  cools.  The  two  models  differ 
in  the  mechanism  that  produces  the  heating.  A  comparison  of 
the  two  models  and  an  evaluation  of  NO  yields  is  given  by  Hill 
(1979b) . 

The  shock  heating  model  is  based  on  the  shock  wave  theory 
of  Lin  (1954).  The  model  assumes  an  infinitely  narrow  line 
source  for  lightning  and  an  instantaneous  release  of  energy 
which  creates  a  sharp  pressure  discontinuity.  Heating  occurs 
as  the  shock  wave  that  is  created  propagates  radially  from 
the  lightning  sources  region.  It  is  this  heating  that  causes 
the  molecular  dissociation  and  eventual  recombination  to  form 


TABLE  2.  PREVIOUS  THEORETICAL  AND  EXPERIMENTAL 

ESTIMATES  OF  NOx  PRODUCTION  BY  LIGHTNING 


Lightning  Model 

References 

NOx/J 

1016 

molec/j 

NOw/f lash 
1025  molec/ 
flash 

Global  NOx 
10^2  gN/yr 

Shock  wave 

4 

2.5 

1.1 

6 

Tuck  1976 

Experimental 

observation 

- 

10 

- 

NOxon  1976 

Hot  channel 

6 

32-40 

Griffinq  1977 

Shock  wave 

3-7 

6-14 

30-40 

Chameides 
et  al.  1977 

Experimental 

spark 

simulation* 

6+1(10)* 

8+2(1350)* 

- 

- 

Chameides 
et  al.  1977 

Shock  wave 

8-17 

16-34 

35-90 

'  Chameides 

1979 

Hot  channel 

80 

6 

Hill  1980 

Review 

- 

3 

Dawsoii  1980 

Experimental 

spark 

simulation* 

5+2 (105)* 

0.5 

1 . 8+0 . 7 

Levine  et  al. 
1981 

Review 

3.7 

0.74 

2.1 

Hameed  et  al. 
1981 

*The  discharge  energy  in  J/m  is  listed  in  parentheses  for  the 
experimental  spark  simulations. 


The  channel  heating  model  is  based  on  theories  of  intense 
localized  sources  of  ohmic  heat  in  the  atmosphere  (Brode  1955; 
Hill  1971;  Plooster  1971).  Lightning  is  simulated  by  a  time- 
dependent  stroke  which  heats  a  channel  of  air  a  few  centimeters 
wide  as  it  propagates.  The  molecular  heating  is  caused  directly 
by  the  lightning  stroke.  At  first  the  extreme  temperatures 
and  high  electron  densities  cause  rapid  molecular  dissociation, 
excitation,  and  ionization  but  as  the  hot  air  expands  and  cools, 
the  ions  and  molecules  recombine. 

2.  Freeze-Out  Temperatures 

Both  mechanisms  result  in  heated  air  that,  presumably,  can 
cause  the  massive  dissociation  and  ionization  of  atmospheric 
species  that  has  been  observed  experimentally  from  emission 
and  absorption  lines  (Salanave  et  al.  1962;  Salanave  1969; 
Wallace  1960;  Orville  1967,  1968).  Early  in  the  lifetime  of 
the  lightning  stroke,  radiative  emission,  collisional  deexci¬ 
tation  and  ionic  recombinations  neutralize  the  plasma.  When 
the  temperature  has  dropped  to  about  4000  K,  neutral  recombina¬ 
tion  reactions  to  form  NO  can  proceed  rapidly.  If  the  cooling 
rates  are  sufficiently  slow,  chemical  equilibrium  can  exist 
for  the  production  and  destruction  of  NO.  When  the  air  has 
cooled  further,  to  the  range  between  2000-3000  K,  the  actual  NO 
concentration  will  not  fall  as  fast  as  the  equilibrium  concen¬ 
tration.  Beyond  this  point  the  NO  formed  at  higher  temperatures 
is  no  longer  efficiently  destroyed  and  it  remains  in  the  air  at 
the  concentration  characteristic  of  the  freeze-out  temperature 
(Gilmore  1974,  1975).  Volume  mixing  ratios  of  a  few  percent  NO 
in  equilibrium  air  can  result. 

The  precise  amount  of  NO  produced  is  very  sensitive  to 
the  freeze-out  temperature  which  in  turn  depends  on  the  cool¬ 
ing  rate  of  the  lightning  channel  gases.  Faster  cooling  rates 
result  in  higher  freeze-out  temperatures.  The  simulations  of 
the  shock  wave  and  hot  channel  models  for  similar  lightning 
initial  conditions  show  that  the  shock  wave  travels  faster  than 


the  heated  gases  in  the  hot  channel.  Therefore,  for  the  shock 
model  calculations  one  would  expect  a  faster  cooling  rate  for 
the  gases  at  the  shock  front  and  a  higher  freeze-out  temperature 
Model  calculations  with  identical  initial  conditions  indicate 
freeze-out  at  3000  K  for  the  hot  channel  model  and  4000  K  for 
the  shock  heating  model  (Hill  et  al.  1980),  confirming  this 
expectation. 

3 .  Model  Comparison 

The  difference  between  the  lightning  models  becomes  most 
evident  when  one  compares  the  volume  of  air  that  has  been 
raised  to  the  freeze-out  temperature.  As  mentioned  previously, 
model  calculations  by  Hill  et  al.  (1980)  predict  that  early 
in  the  lifetime  of  the  lightning  channel,  the  shock  wave  prop¬ 
agates  faster  and  soon  decouples  from  the  hot  gases  in  the 
original  lightning  channel.  This  has  recently  been  confirmed 
by  a  set  of  time-dependent  Schlieren  photographs  of  discharges 
in  air  (Picone  et  al.  1981).  The  experiments  show  that  a  sep- 
ration  of  the  shock  wave  from  the  hot  channel  has  occurred  8 
microseconds  after  the  initiation  of  the  discharge,  and  that  at 
100  microseconds  the  shock  wave  radius  had  quadrupled,  whereas 
the  hot  channel  remained  approximately  the  same  size.  Calcu¬ 
lations  using  both  models  with  identical  initial  conditions 
were  done  to  determine  the  propagation  radius  at  the  freeze-out 
temperature  (Hill  et  al.  1980',.  For  the  hot  channel  model  an 
outer  radius  of  16  cm  was  calculated,  with  freeze^out  occurring 
at  3000  K.  The  shock  wavefront  in  the  shock  heating  model 
had  travelled  1.6  cm  before  freeze-out  at  4000  K.  The  smaller 
volume  of  air  calculated  by  the  shock  heating  model  is  partially 
offset  by  a  greater  air  density.  The  shock  wavefront  envelops 
air  with  ambient  density,  whereas  the  hot  air  density  within 
the  hot  channel  near  the  freeze-out  temperature  is  approximately 
one  tenth  ambient  (i.e.,  300  vs.  3000  K) .  As  a  result  of  these 
considerations  alone,  the  NO  yield  estimates  per  joule  discharge 
energy  for  the  snock  heating  model  are  at  least  a  factor  of 
five  and  perhaps  ten  less  than  the  hot  channel  estimates. 


B.  EXPERIMENTAL  SIMULATIONS  WITH  LABORATORY  DISCHARGES 

Experimental  determinations  of  the  amount  of  N0X  produced 
per  joule  discharge  energy  have  been  studied  in  both  the  high 
arid  low  discharge  energy  regimes.  Chameides  (1979)  pointed  out 
that  it  is  possible  that  results  from  low-energy  discharges  may 
underestimate  the  N0X  production  by  lightning.  An  examination 
of  three  discharge  experiments  over  the  energy  regime  below 
10®  j/m,  as  well  as  various  studies  of  nuclear  explosions,  indi¬ 
cates  an  energy-dependent  trend  for  the  NO  production  efficiency 
which  was  reproduced  by  a  shock  model  calculation  (Chameides 
1979).  The  validity  of  this  trend  needs  to  be  verified  by  an 
experimental  determination  at  higher  energies  other  than  nuclear 
explosions  since  it  is  not  clear  that  the  mechanism  and  time 
scales  of  nuclear  explosions  and  linear  discharges  are  compar¬ 
able.  At  lower  energies,  the  experimental  uncertainty  is  large 
enough  to  make  conclusions  about  a  trend  uncertain. 

1 .  Low-Energy  Discharges 

Lightning  parameters  have  been  studied  with  a  variety  of 
experimental  techniques  in  an  extensive  set  of  experiments  on 
the  Westinghouse  4-meter  spark  chamber  (Xrider  et  al.  1968; 
Dawson  et  al.  1968;  Uman  et  al.  1968,  1970).  These  experiments 
have  been  discussed  in  detail  in  a  recent  review  (Dawson  1980) 
and  will  not  be  repeated  here.  Experimental  spark  simulations 
by  Chameides  et  al.  (1977)  with  low  and  medium  energy  deposi¬ 
tion  resulted  in  similar  N0X/J  yields  within  the  limits  of 
uncertainty.  Table  2  shows  that  these  results  are  also  similar 
to  theoretical  values  calculated  with  the  shock  model  approach. 

2 .  High-Energy  Discharges 

An  experimental  laboratory  discharge  capable  of  generating 
between  10®  and  10®  J/m  discharge  energy  was  used  by  Levine  et 
al.  (1981)  to  simulate  lightning.  Their  simultaneous  measure¬ 
ments  of  NOx  =  NO  +  N02»  NO  and  O3  showed  that  initially  only 
NO  was  produced  by  the  discharge.  Chemical  conversion  to  NO2 
by  NO  +  03-*>N02  +  O2  which  follows  rapidly  in  air  muse  occur 


with  ambient  O3  since  no  enhancement  in  O3  concentration  was 
observed  after  the  discharge.  The  value  obtained  for  the  NOx 
production  per  joule  discharge  energy  at  high  energies  is  in 
line  with  previous  experimental  results  at  lower  energies. 

C.  PRODUCTION  BY  LIGHTNING:  NOx/FLASH 
1 .  Previous  Estimates 

Inherent  in  most  estimates  of  the  amount  of  N0X  produced  in 
a  lightning  stroke  is  the  energy  of  a  typical  lightning  stroke 
in  J/m,  the  length  of  an  average  lightning  stroke,  and  the  addi¬ 
tional  N0X  produced  by  subsequent  strokes  within  a  flash.  These 
lightning  parameters  (listed  in  Table  4  on  page  17),  rather  than 
the  chemical  production  factors,  are  the  largest  sources  of 
variability  in  the  previous  estimates  of  N0X  production. 

Until  very  recently,  a  value  of  10^  j/m  has  been  adopted 
for  the  energy  deposited  by  a  lightning  stroke,  based  on  the 
recommendations  of  Chalmers  (1967)  and  Uman  (1969).  A  recent 
review  by  Hill  (1979a)  extensively  discussed  and  reevaluated 
previous  electrical,  optical,  and  acoustic  measurements  of 
energy  dissipation.  The  new  recommendation  of  104  J/m,  which 
is  an  order  of  magnitude  less,  causes  the  NOx  production  per 
stroke  estimated  by  Hill  et  al.  (1980)  to  also  be  an  order  of 
magnitude  lower. 

Values  suggested  for  the  average  stroke  length  vary  be¬ 
tween  5  and  10  km.  Lightning  channel  lengths  are  typically 
longer  than  5  km  for  tropical  thunderstorms  where  most  light¬ 
ning  activity  occurs. 

Another  parameter  that  needs  to  be  estimated,  but  for 
which  no  direct  experimental  evidence  exists,  is  the  amount  of 
additional  NOx  produced  by  subsequent  lightning  strokes  beyond 
the  first  one.  Chameides  et  al.  (1977)  and  Chameides  (1979) 
consider  an  average  of  four  strokes  per  flash  and  count  each 
one  as  a  separate  event.  Most  probably  this  will  overestimate 
the  NOx  produced.  The  quantity  that  needs  to  be  considered  is 
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the  increase  of  heated  air  in  the  lightning  channel  caused  by 
the  subsequent  strokes  within  the  flash.  Hill  et  al.  (1980) 
assumes  that  the  amount  of  energy  deposited  to  the  lightning 
channel  is  proportional  to  the  charge  deposited  by  each  in¬ 
dividual  stroke  and  that  the  first  stroke  is  about  five  times 
as  energetic  as  the  subsequent  strokes.  For  a  three-stroke 
flash.  Hill  et  al.  (1980)  obtains  heated  air  for  the  entire 
flash  that  is  1.5  +_  0.5  times  the  amount  for  a  single  stroke. 
Dawson  (1980)  favors  even  a  lower  estimate.  He  argues  that 
the  first  stroke  is  by  far  the  most  energetic  and  that  all 
strokes  within  a  flash  occur  along  the  identical  channel  with 
insufficient  time  between  strokes  for  diffusive  losses  of  NOx. 
His  review  recommends  counting  first  strokes  only.  Direct 
experimental  verification  is  needed;  however,  some  insight  may 
be  obtained  from  a  set  of  low-energy  discharge  experiments  by 
Chameides  et  al.  (1977).  The  rate  of  sparking  through  air  was 
varied  between  8  and  53  per  minute.  The  concentration  of  the 
NOx  produced  was  proportional  to  the  total  energy  deposited; 
a  separate  dependence  on  the  sparking  rate  was  not  observed. 
This  result  tentatively  confirms  the  lower  estimates  for  addi¬ 
tional  NOx  produced  within  multiple-stroke  flashes  but  it 
needs  to  be  verified  at  higher  energies  and  at  higher  sparking 
frequencies  to  properly  simulate  the  multiple-stroke  case  for 
lightning. 

2.  In-Situ  NOx  Measurements 

The  only  direct  measurements  of  gas  phase  NCX  produced  by 
lightning  are  the  spectroscopic  determinations  by  Noxon  (1976; 
1978).  Tropospheric  measurements  were  obtained  from  the  ground 
by  observing  spectra  of  the  overhead  sky  during  a  lightning 
storm  with  specific  attention  to  the  wavelength  region  near 
450  nm  that  primarily  contains  NO2  absorption  features.  Ambi¬ 
ent  NO2  concentrations  were  typically  less  than  0.1  ppb.  in 
the  first  experiment  the  large  enhancement  of  over  500  percent 
observed  in  the  early  part  of  the  storm  was  attributed  to 
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lightning  (Noxon  1976).  The  parcel  of  air  with  the  high  NO2 
concentration  was  assumed  to  be  travelling  with  a  storm  cell 
in  which  the  lightning  flash  rate  of  5  per  minute  persisted 
for  one  hour.  A  N02  production  rate  of  1026  N02/flash, 
certain  to  within  one  order  of  magnitude,  was  derived  from  the 
overhead  column. 

It  was  not  possible  to  determine  how  much  of  the  measured 
NO2  was  in  the  clouds  overhead  where  the  absorption  signal 
could  be  greatly  enhanced  by  multiple  scattering.  A  second  set 
of  experiments  in  which  stormy  areas  were  observed  horizontally 
below  the  level  of  the  clouds  resolved  this  difficulty  (Noxon 
1978).  Table  3  shows  a  data  summary  for  experiments  under 
various  lightning  conditions.  The  results  confirm  that  1026 
N02  molecules  are  created  per  flash  beneath  the  clouds.  Obser¬ 
vations  of  the  sky  where  only  corona  discharges  were  occuring 
showed  NO2  concentrations  comparable  to  the  typical  ambient 
values,  indicating  that  they  were  a  minimal  source  of  NO2  during 
thunderstorms.  An  in-situ  measurement  such  as  this  avoids  many 
of  the  large  uncertainties  associated  with  lightning  discussed 
in  the  previous  section  as  well  as  some  uncertainties  in  the 
rainout  process  discussed  in  Chapter  IV.  Noxon's  (1976;1978) 
value  of  1026  N02  molecules  produced  per  flash  from  cloud-to- 
ground  will  be  used  for  further  calculations  in  this  work.  In 
order  to  properly  evaluate  this  recommendation,  a  more  detailed 
discussion  of  the  experimental  uncertainties  follows. 

Although  NO  is  initially  created  in  a  lightning  discharge, 
Noxon's  measurement  of  NO2  approximates  the  NOx  produced,  since 
under  dark  sky  conditions  prevalent  during  a  thunderstorm  most 
of  the  NO  is  rapidly  converted  to  NO2  by  reaction  with  ozone. 

In  full  daylight  the  NO2  is  photolyzed  to  NO;  a  photochemical 
steady  state  exists  and  NO  and  NO2  rapidly  cycle  among  them¬ 
selves.  Immediately  after  a  sudden  decrease  in  solar  intensity, 
such  as  in  an  eclipse,  the  cycling  ceases  and  NO2  builds  up  on 
the  order  of  a  few  minutes  in  the  troposphere.  The  rate  at 


which  this  build-up  occurs  is  limited  by  the  ambient  ozone  con¬ 
centration  and  the  rate  of  mixing  of  the  lightning  channel  air 
with  the  ambient  surroundings.  The  most  serious  uncertainty 
inherent  in  Noxon's  measurements  is  the  possibility  that  the  NO2 
may  have  been  entrained  from  the  stratosphere  where  higher  mix¬ 
ing  ratios  of  NO2  exist  as  the  result  of  the  strong  updrafts 
that  occur  during  thunderstorms.  In  order  to  address  some  of 
these  questions,  in-situ  NO2  measurements  both  below  and  within 
the  clouds  are  needed. 


Table  3.  N02  PRODUCED  RY  LIGHTNING 


Date 

Conditions  N02 

Abundance 

PPb 

Fritz  Peak 

July  1975 

zenith 

cloud,  4  flash/min 

100 

July  1976 

zenith 

cloud,  1  flash/min 

25, <3* 

July  10, 

1976 

zenith 

cloud,  1  flash/min 

20, <3* 

July  18, 

1976 

at  4-km  distant  sunlit  peak 

2 

July  21, 

1976 

zenith 

cloud  after  3  strokes 

1 0  ,< 2t 

July  23, 

1976 

sunlit 

peak,  1  flash/5  min 

3 

July  31, 

1976 

sunlit 

peak,  1  flash/min 

10 

Langmuir  Laboratory 

Aug.  19, 

1976 

sunlit 

horizon  sky,  3  flash/min 

50 

Aug .  21, 

1976 

sunlit 

horizon  sky,  3  flash/min 

30 

Aug.  23, 

1976 

sunlit 

horizon  sky,  corona,  no  lightning 

<0.3 

Aug.  26, 

1976 

sunlit 

horizon  sky,  corona,  no  lightning 

<0.1 

*  1  hour 

*  1  hour 

later, 
earlier . 

Source: 

Noxon 

(1978) 
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D.  PREVIOUS  ESTIMATES  FOR  GLOBAL  NOx 

Table  4  lists  the  assumptions  that  have  been  made  by  pr^- 
l  vious  workers  about  lightning  parameters  to  obtain  the  global 

NOx  production  rates  given  in  Table  ?.  The  accumulation  of 
uncertainties  can  readily  be  seen  in  the  range  of  various  param¬ 
eters,  especially  when  they  are  combined  to  give  the  global  aver 
|  age  of  the  total  energy  dissipated  in  J/cm^s. 

The  annual  number  of  lightninq  flashes  to  ground  that  occur 
worldwide  is  generally  given  as  100/s  and  is  the  value  adopted 
by  all  of  the  works  listed  in  Table  4.  Although  this  estimate 
is  based  on  weak  historical  evidence  of  visual  counting,  more 

i 

recent  measurements  do  not  seriously  contradict  it.  The  review 
by  Orville  and  Spencer  (1979)  indicates  that  the  values  show  a 
seasonal  dependence  and  that  the  average  appears  to  be  less  by 
about  a  factor  of  two.  Comparisons  among  the  different  estimate 
for  NOx  production  must  be  done  carefully  since  not  all  the 
estimates  include  a  contribution  from  intracloud  discharges. 

The  globally  integrated  rate  of  intracloud  discharges  is  gener- 
|  ally  quoted  as  300-400  flashes/s  (Uman  1969;  Chalmers  1967). 

It  must  be  recognized  that  experimental  data  about  intracloud 
discharges  is  generally  obtained  from  indirect  means  such  as 
surface  electric  field,  electric  field  charge,  and  thunder 
measurements;  therefore  the  experimental  uncertainties  are  quite 
large.  Also  the  body  of  information  about  intracloud  discharges 
is  smaller  and  more  recent  than  the  more  easily  studied  ground- 
stroke  phenomenon.  More  about  the  lightning  parameters  and 
their  uncertainties  will  be  said  in  Chapter  3,  where  the  global 
distribution  data  used  in  this  study  are  discussed. 
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TABLE  4.  SUMMARY  OF  LIGHTNING  FLASH  DATA  USED  FOR  GLOBAL 
NOx  ESTIMATES 


Energy 

Deposited 

J/m 

Stroke 

Length 

km 

Strokes 

per 

Flash 

Flashes 

per 

Second 

Total  Energy 
Dissipated 
j/cm2s 

References 

105 

6.25 

1 

100  CG  a 
400  IC 

6.2xl0-8 

Tuck  1976 

105 

5 

4 

100  CG  b 
300  IC 

1.6xl0~7 

Chameides 
et  al.  1977 

105 

5 

4 

100  CG  b 
300  IC 

1.6x10-7 

Chameides 

1979 

104 

5 

1.5 

100  CG  c 

1.5xl0-9 

Hill  1980 

104 

10 

1.0 

100  CG  a 
400  IC 

l.OxlO-8 

Dawson  1980 

- 

- 

- 

- 

10“8 

Levine 
et  al.  1981 

CG  =  Cloud-to-ground.  IC  =  Intracloud. 

a.  Uman  1969 

b.  Chalmers  1967 

c.  Turman  1978;  Orville  and  Spencer  1979. 


III.  GLOBAL  DISTRIBUTION  OF  LIGHTNING  AND  NOx 


A.  GLOBAL  LIGHTNING 

So  far  the  approach  taken  has  been  to  adopt  Noxon's  (1976) 
in-situ  measurement  for  NOx/flash.  Lightning  count  data  and 
global  lightning  distributions  are  required  next.  Most  recent 
global  distribution  data  are  obtained  from  satellites  which 
measure  either  the  direct  flash  rate  with  visible  or  near  IR 
sensors  or  an  indirect  quantity  which  can  be  related  to  flash 
rate,  as  is  done  with  measurements  of  radio  frequency  emission. 
The  detection  threshold  of  the  optical  counting  sensors  needs 
to  be  high  to  eliminate  as  much  noise  as  possible,  so  that 
only  the  most  energetic  lightning  flashes  are  recorded.  In  both 
types  of  measurements  a  normalization  to  an  absolute  flash  rate 
is  necessary. 

1 .  Local  Flash  Rate 

The  recent  review  by  Orville  and  Spencer  (1979)  discusses 
a  set  of  electric  field  lightning  experiments  over  Jacksonville, 
Florida  on  July  22,  1974  for  which  DMSP  satellite  data  was  re¬ 
corded  20  minutes  later.  The  local  flashing  rate  obtained 
on  the  qround  of  approximately  1  s“l  compares  reasonably  well 
with  the  local  rate  of  0.4  s”l  estimated  from  the  satellite  if 
one  assumes  that  the  characteristics  of  the  storm  did  not 
change  much  in  the  intervening  time. 

Electric  field  measurements  by  Livingston  and  Krider  (1978) 
have  been  used  to  normalize  the  most  recently  reported  DMSP 
satellite  data  by  Turman  and  Edgar  (1980).  The  ground  measure¬ 
ments  from  June  1  to  July  15,  1575  were  combined  with  the 
ground  measurements  of  Jacobson  and  Krider  (1976)  between  June 
15  and  July  30,  1974.  In  a  300-km2  experimental  area  at  the 
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Kennedy  Space  Center,  about  38  percent  of  the  lightning  flashes 
were  to  ground  and  the  total  lightning  density  to  ground  was 
about  6  f lashes/km2/month. 

The  use  of  local  flash  rate  data  to  normalize  satellite 
data  brings  to  mind  several  difficulties.  If  the  identical 
storm  is  not  observed  by  both  systems,  sufficient  data  averag¬ 
ing  must  be  done  to  incorporate  the  natural  variations.  The 
detection  systems  used  are  generally  not  the  same;  for  the  two 
examples  discussed  here,  electric  field  measurements  were  done 
on  the  ground  and  optical  sensors  were  used  from  the  satellite. 
The  detection  efficiencies  may  vary  for  ground  flashes  and 
most  probably  do  vary  for  cloud- to-ground  vs.  intracloud 
flashes. 

2 .  Intracloud  vs  Cloud-to-Ground  Lightning 

An  experimental  determination  of  the  ratio  of  the  number 
of  intracloud  discharges  to  the  cloud-to-ground  strokes, 

Nc/Ng,  requires  counting  both  lightning  events  with  equal 
efficiency.  In  practice,  this  is  difficult  to  do  since  the 
ease  of  counting  cloud-to-ground  strokes  tends  to  bias  the 
ratio  to  values  which  are  too  low.  Accurate  lightning  counting 
of  intracloud  discharges  by  optical  means  is  very  difficult. 
First  of  all,  in  many  cases  the  visibility  within  a  cloud  is 
too  poor  to  permit  identification  of  a  lightning  event.  Also, 
precipitation  below  the  clouds  diminishes  the  visibility  even 
further.  At  long  distances,  discharges  within  the  clouds 
often  cannot  be  distinguished  from  ground  flashes. 

Various  indirect  methods  have  been  used  to  distinguish 
between  the  two  forms  of  lightning.  Radio  noise  emitted  from 
lightning  discharges  can  be  used  to  selectively  detect  either 
cloud-to-ground  strokes  or  intracloud  discharges,  depending 
on  the  frequency  region  monitored.  Table  5  shows  that 
emission  from  the  ground  stroke  predominates  at  lower  frequen¬ 
cies  and  at  higher  frequencies,  above  about  1  MH~,  both  emit 
with  equal  intensity.  At  still  higher  frequencies,  greater 


than  about  10  MHz,  Kreielsheimer  and  Lodge-Osborn  (1971)  found 
that  gaps  in  the  radiative  noise  from  the  two  types  of  light¬ 
ning  followed  different  distribution  laws  which  could  be  used 
to  selectively  choose  between  them.  Ground-based  measurements 
can  make  use  of  the  entire  frequency  range,  but  measurements 
from  satellites  must  be  done  in  frequency  regions  that  exceed 
the  critical  frequency  of  the  ionosphere.  Electric  field 
change  measurements  can  also  be  used  to  distinguish  between 
the  two  types  of  lightning  by  the  different  features  present 
in  the  electric  field  change  time  record.  A  more  detailed 
discussion  of  the  experimental  methods  used  to  determine  the 
ratio  Nc/Ng  is  presented  by  Prentice  (1977). 


TABLE  5.  RADIO  FREQUENCY  EMISSION  FROM  CLOUD 
AND  GROUND  STROKES 


Frequency 

Amplitude  Ratio 
(ground  stroke/cloud  stroke) 

3  kHz 

6  kHz 

10  kHz 

20  kHz 

30  kHz 

50  to  100  kHz 

1.5  to  12  MHz 

”  20/1  to  40/1 

10/1  to  20/1 

10/1 

5/1 

2/1  to  3/1 

1/1  to  1.5/1 

1/1 

Source:  Golde  1977,  p.  223. 

Prentice  and  Mackerras  (1977)  have  analyzed  29  data  sets 
from  13  countries  and  obtained  the  following  relationship, 
which  generalizes  the  ratio  Nc/Ng  in  a  typical  thunderstorm, 

n  ,  .  0.4T  1  ,,  . 

0.6  72  -  0 . 98  A  J  (3,1 


Nc  t  a  ic 
Ng  =  (4‘16 


+  2.16  cos  3  A) 


where 


T  £  84,  thunderdays  per  year 

A  <  60°,  latitude  in  Northern  and  Southern  Hemispheres. 


The  authors  recommend  equation  3.1,  but  if  data  for  both  the 
annual  number  of  thunderdays  and  the  latitude  are  unknown,  the 
following  relationships  can  be  used  to  estimate  the  ratio 
Nc/Ng  from  either  the  annual  thunderdays  or  the  latitude. 

SS  -  1.0  +  0.063  T,  (10<  T  <  84)  (3.2) 

Ng  —  — 

Nc 

^  =  4.16  +  2.16  cos  3 A ,  (0  <x<  60°)  (3.3) 

The  ratio  as  a  function  of  latitude  from  equation  3.3  is  shown 
in  Figure  1.  This  ratio  has  the  highest  value  in  the  ttopics, 
where  most  lightning  activity  occurs. 

The  partitioning  between  intracloud  discharges  and  cloud-to- 
ground  lightning  can  be  expressed  in  fractional  form  as  follows: 

c  .  Ng  _  _  _ 1 

9  ”  Ng+NC  Nc/Ng  +  1 

f  Nc  Nc/Ng 

tc  Ng+Nc  Nc/Ng  +  1 

The  amount  of  NOx  formed  by  intracloud  discharges  and 
cloud-to-ground  strokes  can  be  estimated  from  the  energy  re¬ 
leased  by  intracloud  lightning  processes  and  cloud-to-ground 
lightning  processes.  The  energy  of  a  cloud- to-ground  stroke 
can  be  represented  by  E  =  1/2  QU  where  Q  is  the  charge  trans¬ 
ferred  and  U  is  the  electrostatic  potential.  It  is  not  possible 
to  estimate  the  potential  of  a  thundercloud  with  respect  to 
ground  because  the  charge  distribution  within  a  cloud  is  unknown; 
however,  it  is  possible  to  estimate  the  potential  of  the  upper 
end  of  the  leader  and  ground  before  it  reaches  the  earth.  In 
a  recent  review  by  Berger  (1977),  typical  values  were  8  coulombs 
for  the  charge  and  60  V/cm  for  the  gradient  along  the  leader 
«hich  becomes  30  MV  for  a  5-km  stroke.  The  voltage  from  the 
accompanying  corona  shell  was  added,  resulting  in  a  total 
voltage  of  about  50  MV.  Berger  concluded  that  a  typical  cloud- 


(3.4) 

(3.5) 
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FIGURE  1.  Relation  between  observed  values  of  the 
ratio,  Nc/Ng  of  cloud  flashes  to  ground 
flashes  and  latitude,  The  29  data 
points  and  the  empirical  relation 
Nc/Ng  =  4.16  +  2.16  cos  3  \  a.-e  shown. 

A  similar  analysis  can  be  done  for  the  K  change  stage  of 
an  intracloud  discharge  which  is  analogous  to  the  return  stroke. 
The  review  by  Brook  and  Ogawa  (1977)  lists  a  range  of  dissipated 
charge  of  10  to  100  coulombs  with  a  typical  value  of  about  50 
coulombs.  A  typical  discharge  length  was  2  km  and  an  average 
electric  field  change  upon  discharge  was  100  V/m.  If  one 
assumes  a  simple  dipole  interaction  for  the  discharge  process, 
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the  energy  deposited  by  intracloud  discharges  is  at  least  one 
order  of  magnitude  less  than  the  typical  cloud-to-ground  stroke. 
It  should  be  noted  that  the  electric  field  change  measurements 
may  be  indicative  of  charge  neutralization,  but  it  can  also  be 
caused  by  charge  rearrangement  within  the  cloud.  Typical  cur¬ 
rent  estimates  also  support  the  values  estimated  for  the  depos¬ 
ited  energy;  the  current  dissipated  by  a  typical  intracloud 
discharge  was  estimated  by  Brook  and  Ogawa  (1977'  to  be  1400 
amperes  whereas  the  current  of  a  cloud-to-ground  stroke  is 
typically  ten  times  larger.  Both  the  potential  and  current 
arguments  suggest  that  cloud-to-ground  strokes  are  about  ten 
times  more  energetic  than  cloud  discharges. 

Holmes  et  al.  (1971)  reviewed  several  experimental  deter¬ 
minations  of  acoustic  energy  released  by  intracloud  and  cloud- 
to-ground  lightning.  Integrated  thunder  power  spectra  indicated 
that  intracloud  discharges  are  one  third  as  energetic  as  ground 
flashes. 

There  is  a  large  uncertainty  in  all  the  estimates  of  light¬ 
ning  energy  release.  The  deposition  of  energy  into  air  is  gen¬ 
erally  inferred  from  indirect  evidence  and  is  difficult  to  esti¬ 
mate  without  ambiguity  since  the  exact  mechanism  of  energy  dep¬ 
osition  is  not  currrently  understood.  For  example,  the  theory  o 
the  frequency  of  spectrum  of  thunder  is  still  a  matter  of  contro 
versy,  so  that  the  use  of  integrated  power  spectra  to  estimate 
the  energy  release  is  uncertain.  These  problems  are  extensively 
discussed  in  Hill's  (1979)  survey  of  lightning  energy  estimates. 

In  this  report,  cloud-to-ground  strokes  will  be  considered 
to  be  ten  times  more  efficient  as  NOx  producers  than  intracloud 
discharges  since  we  assume  that  cloud-to-ground  strokes  deposit 
ten  times  more  energy  per  stroke  on  the  average.  If  Holmes' 
factor  of  3  were  used  instead,  the  results  would  show  higher 
N0X  values  at  higher  altitudes,  but  the  total  NOx  production 
would  only  increase  about  50  percent,  well  within  the  estimated 
uncertainty. 
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3 .  Global  distribution  of  Lightning 

There  are  a  number  of  measurements  of  the  global  distribu¬ 
tion  of  lightning: 

(1)  Optical  observations  of  lightning  flashes  from  the 
various  DMSP  satellite  experiments  by  Orville  and 
Spencer  (1979)  and  Turman  and  Edgar  (1980)  -  see 
Figure  2. 

(2)  Observations  of  HF  radio  noise  from  the  ISS-b 
satellite  [Kotaki  et  al.  (1981)]  -  see  Figure  3. 

(3)  Observations  of  high  frequency  radio  noise  from 
the  ground  by  Davies  (1965)  -  see  Figure  4.,  also 
Crichlow  et  al.  (1971). 

(4)  Surface  estimates  of  thunderstorm  days  compiled  in 
the  Handbook  of  Geophysics  (1966)  -  see  Figure  5. 

Overall,  these  results  seem  to  be  fairly  consistent,  ex¬ 
cept  for  the  absolute  variation  by  a  factor  of  3  between 
items  1  and  2.  The  strengths  and  weaknesses  of  the  data  sets 
are  listed  in  Table  6. 

The  following  general  characteristics  are  apparent: 

o  Lightning  occurs  predominantly  over  or  near  land,  at 
low  or  middle  "meteorological"  latitudes  [namely,  the 
latitudes  relative  to  the  Inter-Tropical  Convergence 
Zone  (ITCZ),  whose  mean  position  is  at  5-10°N  geo¬ 
graphic  latitude  year-round,  but  with  relatively 
large  seasonal  (Monsoonal)  excursions,  especially 
over  the  Indian  Ocean].  There  are  very  few  thunder¬ 
clouds  at  high  latitudes  because  solar  forcing  is 
very  weak  there,  nor  over  the  open  ocean  where  the 
surface  is  so  uniform  that  few  large  disturbances 
are  likely  to  be  launched  up  to  become  cumulonimbus 
clouds. 

o  A  longitudinal  summation  of  Turman  and  Edgar's  (1980) 

DMSP  data  is  shown  by  the  bar  graphs  in  Figure  6.  It 
clearly  shows  that  at  dusk  most  of  the  lightning  activity 
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LONGITUDE,  deg 

FIGURE  2.  Optical  observation  of  lightning  rate  from  DMSP 
satellite  for  dusk  in  April.  (Turman  and  Edgar 
1980) 
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FIGURE  3.  Lightning  rate  data  from  the  HF  radio  noise  data 
obtained  by  the  ISS-b  satellite  for  spring  in  the 
late  afternoon  and  evening  (Kotaki  et  al.  1981) 


FI CURE  4.  Global  distribution  of  atmospheric  radio  noise 

obtained  from  ground  monitoring  stations.  Spring 
season,  1200-1600  hours  (CCIR  Report  322  1963) 
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FIGURE  5.  Springtime  global  distributions  determined  from 
ground-based  observations  of  thunderstorms  (Hand 
book  of  Geophysics  1960) . 


TABLE  6.  STRENGTHS  AND  WEAKNESSES  OF  THE  PRIMARY  DATA  SETS 


Data  Set 

Strengths 

Weaknesses 

DMSP 

Satellite  provides 
good  global  coverage. 

Data  only  at  dawn 
and  dusk. 

A  flash  is  a  good 
direct  indicator. 

Calibration:  what 

is  flash  threshold? 

ISS-b 

Satellite  provides 
good  global  coverage. 

Calibration:  the  fac¬ 

tor  of  three  between 
DMSP  and  ISS-b  needs 
to  be  reconciled. 

Radio  noise  penetrates 
clouds,  so  in-cloud 
flashes  are  counted. 

possible 
interferences . 

Pulse  analysis  made 
possible  direct 
counting  of  lightning 
events . 

HF  Radio 

Noise 

Long-term  data  base. 

Ground  based:  sparse 
global  coveraqe-- 
16  stations. 

Calibration:  quan¬ 

titative  relation¬ 
ship  between  radio 
noise  intensity  and 
the  lightning  rate 
is  unknown. 

Thunderstorm 

Day 

Long-term  data  base. 

Ground-based,  i.e., 
not  very  uniform 
coverage . 

— 

Calibration:  what 

is  the  relation  be¬ 
tween  thunderstorm 
day  and  lightning 
flash  rate? 
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LATITUDE 

FIGURE  6.  Latitudinal  and  seasonal  dependence  of  the  light¬ 
ning  rate  observed  from  the  DMSP  satellite. 

(Turman  and  Edgar  1980) .  The  curve  shown  for  the 
dusk  data  was  calculated  by  equation  3.6. 


occurs  within  the  latitude  range  +30.  Seasonal  excur¬ 
sions  of  peak  lightning  activity  into  the  midlatitudes 
occurs  during  the  local  summer,  as  one  would  expect. 

•  On  a  diurnal  basis  the  frequency  of  lightning  tends 
to  peak  in  the  late  afternoon. 

A  simple  gaussian  analysis  of  the  seasonal  and  latitudinal 
behavior  is  given  by  the  solid  curve  in  Figure  6  and  the 
resulting  mathematical  expression  for  lightning  flash  rate  in 
flashes  per  second  is, 


L(X ,M) 


l'.Txl'O4 


X0)  /288 


+  (10  +  4A)  e-(X-35)Vl28 


(3.6) 


where 

X  »  16  cos  l 30( M-7) ] 

A  =  cos  (30(M-6) ] 

=  latitude,  include  negative  sign  for  Southern 
Hemisphere 

M  =  month  number,  i.e.,  Jan  =  1,  Feb  =  2,  etc. 

L  =  lightning  rate,  annual  average  in  flashes/s. 

It  shows  a  small  peak  around  35°N  latitude  corresponding 
to  lightning  off  the  eastern  coasts  of  the  Southern  U.S.  and 
Japan  that  persists  during  all  seasons.  During  the  Northern 
Hemisphere  summer  the  larger  peak  is  predominantly  caused  by 
lightning  activity  over  North  America  and  the  Indian  Ocean  area, 
whereas,  during  the  Southern  Hemisphere  summer  the  peak  light¬ 
ning  activity  occurs  in  Brazil  and  the  southern  half  of  Africa. 

Table  7  summarizes  the  latitudinal  and  monthly  behavior 
of  the  lightning  rate  calculated  from  equation  3.6.  It  should 
be  noted  that  equation  3.6  merely  partitions  the  total  light¬ 
ning  rate  among  the  months  of  the  year  and  latitudes  between 
+60°,  so  that  this  treatment  is  valid  for  any  global  lightning 
rate.  A  value  of  300  flashes/s,  including  both  intracloud 
discharges  and  cloud-to-ground  flashes,  ha's  been  assumed. 
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TABLE  7.  LIGHTNING  FLASH  RATE,  LATITUDINAL  AND  SEASONAL  DISTRIBUTION* 


An  average  global  flash  rate  of  300  flashes/second  has  been  assumed.  The  values  listed  have  been 
calculated  for  each  degree  of  latitude,  and  summed  over  10-degree  latitude  bins. 


4 •  Globally  Integrated  Lightning  Rate 

Globally  integrated  rates  of  lightning  observed  by  various 
satellite  experiments  are  summarized  in  Table  8.  Differences 
among  the  experiments  can  be  expected,  since  some  experiments 
were  done  at  dawn  and  dusk  and  others  at  midnight.  Most  strik¬ 
ingly,  the  only  experiment  not  done  with  an  optical  sensor 
records  about  three  times  more  lightning  than  the  previous 
works  (Kotaki  et  al.  1981). 

The  geographical  distribution  for  the  ISS-b  and  the  most 
recent  DMSP  experiment  are  similar  (see  Figures  2  and  3).  It 
is  possiole  that  the  difference  in  the  flashing  rate  results 
from  the  method  used  to  calibrate  the  satellite  data.  The 
following  discussion  summarizes  the  experimental  techniques 
and  calibration  methods  used  to  obtain  an  absolute  value  for 
the  lightning  flash  rate. 

The  high-frequency  radio  noise  observed  by  the  ISS-b  satel¬ 
lite  (Kotaki  et  al.  1981)  which  was  attributed  to  lightning  was 
found  to  coincide  with  the  presence  of  cumulonimbus  clouds 
observed  at  the  same  time  by  the  geostationary  meteorological 
satellite  GMS.  Radio  noise  data  was  collected  at  2.5,  5,  10, 
and  25  MHz.  The  frequencies  were  chosen  to  be  high  enough  to 
be  above  the  critical  frequency  of  the  ionosphere  yet  low  enough 
to  take  advantage  of  the  greater  signal  strengths.  At  the  same 
time,  it  was  necessary  to  avoid  interference  from  radio  trans¬ 
mitters  . 

The  noise  patterns  were  distinctive  and  were  easily  dis¬ 
tinguished  from  cosmic  radio  noise  which  appeared  as  a  constant 
background  level  and  from  solar  radio  noise  originating  from 
solar  bursts  which  appeared  as  a  sinusoidal  signal  synchronized 
with  the  spin  of  the  satellite.  About  30  percent  of  the  time, 
interference  from  transmitters  on  the  ground  was  evident  but 
its  characteristic  random  noise  pattern,  which  persists  over 
long  periods  of  time,  could  be  readily  distinguished  from  the 
intermittent  waveforms  characteristic  of  the  noise  from  lightning 


TABLE  8.  LIGHTNING  FLASHES  COUNTED  BY  SATELLITE 


Times  during  which  interference  was  evident  were  disregarded. 
Digital  counts  were  obtained  when  the  analog  signal  received  by 
the  satellite  detector  exceeded  a  threshold  level  and  had  the 
correct  waveform.  A  lightning  event  was  said  to  occur  if 
digital  counts  were  observed  in  all  the  frequency  channels  which 
exceeded  the  critical  frequency  of  the  ionosphere  simultaneously. 
This  is  a  direct  digital  counting  method  that  does  not  rely 
on  a  calibration  for  the  absolute  lightning  flash  rate. 

The  DMSP  satellite  measurements  by  Turman  and  Edgar  (1980) 
were  done  optically  with  a  silicon  photodiode  with  a  response 
function  between  500  and  1100  nm.  The  sensor  was  operated 
in  a  triggered  mode  where  only  the  high-frequency,  high-intensity 
pulses  were  recorded  in  order  to  discriminate  against  interfering 
light  sources  from  the  earth.  The  result  was;  that  only  the  most 
energetic  lightning  strokes  were  observed,  creating  a  selective 
discrimination  against  the  less  energetic  intracloud  lightning. 

A  calibration  was  necessary  to  obtain  an  absolute  value  for 
the  lightning  rate. 

Based  on  the  results  of  Forrest  (1950),  Turman  and  Edgar 
(1980)  assumed  that  the  mean  diurnal  flash  rate  is  equivalent 
to  the  instantaneous  rate  at  dusk.  Normalization  for  the 
DMSP  satellite  data  was  obtained  by  comparison  with  the 
ground-based  electric  field  measurements  of  Livingston  and 
Krider  (1978;  also  see  Section  III.A.l)  and  with  various  other 
local  determinations  of  absolute  flash  frequency  to  ground 
(Wanasalyja  1979;  Prentice  1977).  The  total  lightning  flash 
rate  was  inferred  from  Livingston  and  Krider's  (1978)  measure¬ 
ments  for  Eastern  Florida  in  which  about  38  percent  of  the 
lightning  flashes  were  striking  ground.  An  examination  of 
Figure  1  shows  that  in  tropical  regions,  where  most  lightning 
occurs,  the  ratio  Nc/Ng  is  generally  higher. 

If  one  uses  an  intracloud  to  cloud-to-ground  ratio  rang¬ 
ing  between  4  and  6  to  renormalize  the  DMSP  data  for  latitudes 
between  30°N  and  30°S,  about  250  +  50  flashes/s  tesult.  This 


then  agrees  more  closely  with  the  ISS-b  satellite  measurement 
of  about  280  flashes/s. .  In  this  work  a  rounded  value  of  300 
flashes/s  was  adopted  for  the  globally  integrated  lightning 
rate. 

B.  INJECTION  HEIGHT  OF  N0X 

The  discussion  of  injection  height  can  be  conceptually 
considered  in  two  domains,  first  the  altitude  of  the  NOx  pro¬ 
duction  and  second  the  altitude  at  which  the  NOx  is  eventually 
deposited,  which  may  be  different  due  to  the  violent  updrafts 
and  downdrrfts  associated  with  thunderstorms.  Updrafts  of 
60  m/s  in  the  center  of  the  thunderstorm  cloud  surrounded  by 
downdrafts  of  35  m/s  have  been  observed.  Also,  clouds  rise 
buoyantly,  causing  a  net  transport  of  air. 

For  cloud-to-ground  flashes,  the  initial  production  alti¬ 
tude  is  between  the  lower  portion  of  the  cloud  and  the  ground. 
The  thickness  and  altitude  of  thunderclouds  vary  greatly,  but 
lightning  rarely  occurs  in  clouds  less  than  3  km  thick  and 
most  lightning  activity  occurs  in  large  convective  storms 
which  can  extend  up  to  20  km  in  altitude.  Some  mean  estimates 
of  production  height  are  listed  in  Table  9,  but  this  does  not 
directly  answer  the  question  of  effective  injection  height, 
which  depends  on  the  cloud  stabilization  height  rather  than  on 
the  height  of  the  lightning  discharges. 

For  strong  midlatitude  thunderstorms  it  is  well  known  that 
the  anvil  top  of  the  cumulonimbus  cloud  penetrates  into  the 
stratosphere  and  the  extent  of  penetration  has  been  related  to 
the  kinetic  energy  associated  with  updrafts  within  a  thunder¬ 
cloud  (vonnegut  and  Moore  1958).  Penetrations  as  much  as  6  km 
have  been  reported  for  large  thunderstorms  over  Oklahoma  (Roach 
1967).  The  consequences  of  cloud-top  penetration  is  still  a 
matter  of  controversy.  For  example,  it  is  not  clear  whether 
cloud-top  penetration  ultimately  results  in  injection  or  re¬ 
moval  of  chemical  species  such  as  water  vapor  or  ozone  from 
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TABLE  9.  SOME  REPRESENTATIVE  HEIGHTS  OF  CLOUD  DISCHARGES 


He ight 
Range 

Location 

Reference 

5.8  to  4.7 

Albuquerque,  NM 

Workman  and  Holzer  (1942) 

5.8  to  5.2 

Albuquerque,  NM 

Workman  et  al.  (1942) 

5.5  to  5.1 

Socorro,  NM 

Reynolds  and  Neill  (1955) 

>8 

South  Africa 

Malan  (1956) 

10  to  6 

- 

Hatakeyama  (1958) 

6 

- 

Tamura  et  al.  (1958) 

7  -  lla 

3  -  6b 

- 

Takagi  (1961) 

8.2  to  4.8 

Tropics 

Wang  (1963) 

6.0  to  4.0 

Socorro,  NM 

Ogawa  and  Brook  (1964) 

4  to  1 2a 

1  to  8b 

Brisbane,  Australia 

Mackerras  (1968) 

8 

- 

Nakano  (1973) 

5  +  ! 

Socorro,  NM 

Tucson,  AZ 

Teer  and  Few  (1974) 

Roswell,  NM 

Houston,  TX 

ahigh  altitude  intracloud  discharges 

^typical  intracloud  discharges 

Source:  Adapted  from  Brook  and  Ogawa  (1977). 

the  stratosphere  (Barrett  et  al.  1972;  Johnston  and  Solomon 
1979).  Stratospheric  water  vapor  concentrations  near  the 
tropopause  that  are  not  entirely  consistent  with  those  one 
would  expect  from  the  Hadley  cell  global  circulation  model 
also  indicate  that  the  understanding  of  the  mechanism  of  water 
vapor  transfer  across  the  tropopause  is  currently  incomplete 
(Ellsaesser  et  al.  1980).  The  fate  of  NO2  near  the  tropopause 
is  linked  to  the  general  mass  transfer  behavior  but  not  neces¬ 
sarily  to  the  behavior  of  water  vapor,  which  may  be  unique  due 
to  the  existence  of  both  liquid  and  vapor  phases  under  atmos¬ 
pheric  conditions.  If  the  net  result  of  a  large  thunderstorm 
is  injection  of  tropospheric  gases  into  the  stratosphere, 
lightning  would  be  a  source  term  for  the  stratospheric  NOx 
budget  and,  conversely,  if  the  net  result  is  freeze-out  and 
entrainment  from  the  stratosphere,  lightning  would  represent 
a  stratospheric  NOx  sink. 

Even  less  information  is  available  for  tropical  thunder¬ 
storms.  Cloud  tops  as  high  as  20  km  have  been  reported  in  both 
the  tropics  and  at  midlatitudes,  yet  these  values  are  excep¬ 
tional.  In  this  work  we  shall  assume  that  updrafts  will  carry 
chemical  species  to  the  cloud  tops  and  that  the  cloud  tops  fol¬ 
low  the  height  of  the  local  tropopause.  The  following  injection 
heights  have  been  assumed: 

Tropical  injections,  between  +30°  latitude 
intracloud  10  to  15  km 

cloud-to-ground  0  to  10  km 

Midlatitude  injections,  between  +(30  to  60)°  latitude 
intracloud  7  to  12  km 

cloud-to-ground  0  to  7  km 

The  vertical  distribution  within  the  two  discharge  regions 
for  each  latitude  interval  depends  on  the  initial  vertical  dis¬ 
tribution  formed  by  the  lightning  process  and  ultimately  on  the 
turbulent  mixing  that  occurs  during  the  thunderstorm. 
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Equilibrium  calculations  at  elevated  atmospheric  tempera¬ 
tures  show  an  equivalent  ( 1  to  2  percent)  mole  fraction  of  NO 
below  3000  K  for  all  tropospheric  c.ir  densities  (Gilmore  1974; 

Zinn  and  Sutherland  1975).  This  indicates  a  NO  yield  propor¬ 
tional  to  the  local  density  or  P/T.  This  physical  picture  is 
consistent  with  a  density-weighted  vertical  distribution  for 
the  NOx  produced  by  lightning,  assuming  no  other  sources  or 
sinks  on  this  time  scale. 

There  is  little  experimental  verification  of  this  point; 
to  date  all  laboratory  simulations  of  lightning-type  discharges 
in  air  have  been  done  at  one  atmosphere  pressure.  Previous 
photographic  evidence  of  the  lightning  return  stroke  does  indi¬ 
cate  somewhat  more  luminosity  near  the  ground,  but  the  effect 
is  not  very  marked  and  is  difficult  to  quantify  since  some  inten¬ 
sity  reduction  with  height  can  be  expected  from  the  branching 
evident  in  the  return  stroke.  Very  recent  return  stroke  lumi¬ 
nosity  measurements  by  Master  et  al.  (1981)  indicate  that  the 
current  dissipation  is  greatest  near  the  ground  and  drops  off 
with  height  more  rapidly  than  was  previously  suspected  (Uman, 
private  communication).  The  analysis  of  this  data  is  still  in 
progress,  thus  the  effect  is  not  yet  quantified,  but  it  is  clear 
that  a  density  weighting  is  conservative. 

In  this  work,  density-weighted  injection  rates  were  proposed 
for  direct  use  in  a  2-D  atmospheric  model.  The  values  listed 
in  Table  S-2  in  the  Executive  Summary  were  calculated  from  the 
column  rates  in  Table  10  using  atmospheric  densities  from  the 
1976  U.S.  Standard  Atmosphere,  included  on  the  right  side  of 
Table  S-2. 

Additional  measurements  of  cloud  tops  are  necessary  to  re¬ 
duce  the  uncertainty  in  injection  height  estimates  and  to  elucidate 
the  mechanism  of  tropopause  penetration.  A  possible  way  to 
establish  the  variation  of  cloud  tops  with  latitude  is  to  use 
stereoscopic  viewing  of  the  overlap  region  between  the  two  U.S. 

GOES  geostationary  satellites  (W.  Shenk,  private  communication). 
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C.  RECOMMENDED  TREATMENT  OF  LIGHTNING  AND  NOx  DISTRIBUTIONS 
How  should  the  global  distr ibu^ ion  data  be  nodelled?  For 
convenience,  a  summary  of  the  recommended  values  needed  to 
calculate  the  global  distribution  is  included  in  Table  10. 

The  ratio  Nc/Ng  has  been  calculated  from  equation  3.3  and  the 
fractions  of  lightning  to  ground  and  within  the  clouds,  fg 
and  fc  were  calculated  for  10-degree  latitude  intervals  using 
equations  3.4  and  3.5.  N0X  production  by  a  intracloud  discharge 
has  been  assumed  to  be  one  tenth  as  efficient  as  a  cloud-to- 
ground  stroke  based  on  the  discharge  energy  considerations 
discussed  in  Section  III. A. 2. 

The  total  NOx  injection  rate  in  grams  of  nitrogen  per  year 
can  be  calculated  for  each  latitude  and  month  as  follows: 

R(A,M)C  =  fc  L(XrM)  0.1C  (3.7) 

R(A,M)g  =  fg  L(  A  ,  M)  C  (3.8) 

where 

c  =  (lO™  N02^  (  mole  ^  .  /l4  gN\  .  fz.2  x  107s^ 

\  flash  /  \6.02  x  10*^  N02/  \mole  /  \  yr  / 

and  L ( A , M )  is  the  lightning  rate  calculated  from  equation  3.6. 

The  annual  average  injection  rate  is  1.85  Tg  N/yr  for  intra¬ 
cloud  discharges  and  3.82  Tg  N/yr  for  cloud-to-ground  strokes. 

A  depth  of  5  km  was  chosen  for  the  region  where  the  NOx  had 
been  deposited  by  intracloud  discharges  and  the  maximum  altitude 
of  this  region  was  chosen  to  model  the  annual  average  tropopause 
height.  Figure  7  shows  the  vertical  and  latitudinal  distribution 
of  NOx  production  for  the  intracloud  and  cloud-to-ground 
regions.  An  atmospheric  density-weighted  partitioning  with 
altitude  for  each  region  is  suggested.  Values  for  a  2-D  N0X 
distribution  are  listed  in  Table  S-2  of  the  Executive  Summary. 


For  a  3-D  distribution  the  reader  is  referred  directly  to 
the  DMSP  satellite  maps  which  contain  latitudinal  and  longi¬ 
tudinal  lightning  rate  distributions. 


XV.  RAINOUT  OF  NOx  PRODUCED  BY  LIGHTNING 

Historically,  the  search  for  experimental  evidence  of  N0X 
production  by  a  lightning  flash  has  been  approached  by  looking 
for  enhancement  of  the  local  HNO3  concentration  in  rainwater 
collected  on  the  ground  which  was  correlated  with  the  number 
of  flashes  or  the  intensity  of  flashes.  To  date,  no  such  en¬ 
hancement  has  been  documented  (e.g.,  Viemeister  1960;  Visser 
1961;  Reiter  1970).  A  positive  correlation  has  been  found 
between  the  N03“  concentration  dissolved  in  rainwater  and 
both  in-cloud  turbulence  and  the  time  integral  of  electrical 
field  strength  (Reiter  1970  and  references  cited  therein). 

The  NOx  formation  was  attributed  to  corona  discharges  during 
the  thunderstorm.  Experiments  of  this  type  are  difficult  to 
extrapolate  to  gas  phase  NOx  concentrations  in  the  troposphere 
since  dissolved  N03“  is  the  result  of  homogeneous  and  hetero¬ 
geneous  chemical  reactions  that  are  not  well  understood, 
especially  under  thunderstorm  conditions. 

Loss  mechanisms  for  NO2  formed  by  a  lightning  flash  during 
a  thunderstorm  include  the  following: 

•  Gas  phase  reactions  among  N0X  species  that  produce 
HNO3  or  HNO2  during  the  initial  stages  after  the 
lightning  flash  when  high  temperatures  and  high 
local  NOx  concentrations  are  present. 

•  The  gas  phase  reactions  to  form  nitric  and  nitrous 
acids 

N02  +  OH  +  M - ►  HNO3  +  M  (4.1) 

NO  +  OH  +  M - -HN02  +  M.  (4.2) 

This  proceeds  best  in  the  presence  of  sunlight  when 
higher  OH  concentrations  are  present. 


•  The  gas  phase  reaction  sequence, 


(4.3) 

(4.4) 

(4.5) 


NO2  +  O3  —  ■  NO3  +  O2 
NO3  +  N02  +  M  ~ — +  M 
N205  +  H20 - -2  HNO3. 


This  sequence  proceeds  best  in  the  absence  of  sunlight, 
when  photolysis  of  N02»  NO3,  and  N2O5  is  slow. 

•  The  dry  deposition  of  gas  or  particles. 

•  Heterogeneous  reaction  on  particles  or  water  droplets. 
The  lifetime  of  NO2  with  respect  to  each  of  these  processes  is 
listed  separately  in  Table  11.  Since  the  mean  lifetime  of  a 
cumulonimbus  cloud  is  about  30  minutes,  it  is  evident  that  none 
of  these  mechanisms  is  an  important  loss  term. 

TABLE  11.  APPROXIMATE  LIFETIME  FOR  INDIVIDUAL  N02  LOSS  MECHANISMS 


Loss  Meci.nnism 

N O2  Lifetime 

Condition  for  \tolidity 

*1 

N02  +  OH  +  M - -HNO3  +  M 

3  -  17  hrs 

10  hrs 

With  daytime  ambient 
[OH)  -  (1-6)  x  106  air 3 

In  absence  of  sunlight 

N02  ♦  O3  ■*  NOj  +  02 

NOj  +  NOj - *■  N2O5 

k4 

N2°5  +  HgO-  «  2  HNO3 

50  -  500  hrs 

In  absence  of  sunlight 

[H2OJ  -  3  x  1016  an-3 

Dry  deposition 

N02— —  ground 

1-10  days 

Probably  closer  to  1  day 
within  the  planetary 
boundary  layer. 

Heterogeneous  reaction  on 
particles  or  water  droplets 

'* 

•Dissolution  of  gaseous  f©2  81x5  into  water  is  minimal.  The  Henry's  Law 
Coefficients  are  1.2  x  10-2  and  6  x  1(T3  H/atm,  respectively.  Heterogeneous 
reaction  rates  on  particles  depend  on  the  local  aerosol  concentration  levels 
which  vary  widely. 

*1  [M]  ■  1.6  x  10"11  cm3  rolec'i 

k. <10-20  an3  molec-3 


A .  N0X  LOSS 

1.  NOx  Loss:  HNO-t  Formation  in  High  Local  NOx  Air 

Immediately  after  the  lightning  flash  it  is  possible  for 
chemical  conversion  of  N0X  to  HNO3  to  occur  readily  because 
of  the  locally  high  NO  concentrations.  When  the  air  in  the 
lightning  channel  has  reached  the  NO  freeze-out  temperature, 
the  local  NO  mole  fraction  can  be  as  high  as  1  to  5  percent. 
Under  such  conditions,  homogeneous  and  heterogeneous  reactions 
such  as 


2  NO  +  02 - 2  N02  (4.6) 

NO  +  N02  +  H20— «*2  HN02  (4.7) 

NO  +  NO3 - 2  N02  (4.8) 

as  well  as  reaction  sequence  4.3  to  4.5  can  be  extremely  fast. 
As  the  lightning  channel  continues  to  expand  and  mix  with  the 
surrounding  air,  the  temperature  approaches  ambient  and  the 
N0X  concentration  approaches  the  parts  per  billion  range. 

At  complete  dilution  the  rate  of  reaction  among  the  NOx 
reactants  becomes  extremely  slow. 

Another  mechanism  for  nitrous  and  nitric  acid  formation 
has  been  proposed  by  the  time-and-temperature  dependent,  chemi¬ 
cal  kinetic,  hot-channel— model  calculations  of  Hill  and  Rinker 
(1981)  which  focused  on  the  production  of  HNO3  and  HN02  within 
the  heated  air  produced  by  lightning.  The  24  chemical  species 
created  by  lightning  from  N2,  02 ,  H20,  and  C02  were  assumed  to 
be  in  thermodynamic  equilibrium  by  the  time  the  hbt  channel  of 
air  had  expanded  and  cooled  to  3000  K.  The  calculated  time  and 
temperature  dependence  of  HN02  and  HNO3  formation  as  the  gases 
cooled,  expanded,  and  mixed  with  the  surrounding  air  showed  an 
increased  production  that  was  not  simply  fixation  of  products 
formed  at  higher  temperatures.  HN02  and  HNO3  formation  by  re¬ 
actions  4.1  and  4.2  represent  loss  terms  for  the  NO,  N02  and  OH 
concentrations  which  are  maintained  by  thermal  equilibrium,  and 
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thus  a  net  production  of  nitrogen  species  occurs  even  below  the 
temperature  of  NO  concentration  freeze-out.  The  amount  of 
nitrate  formation,  1.6  x  10*5  NO^-  per  cm3  of  heated  lightning 
channel  air  at  3000  K,  was  ccrsistent  with  the  0.2  ppm  N 03“  ob¬ 
served  in  rainfall  at  sea  by  Gambell  and  Fisher  (1964),  assuming 
that  all  the  HNO3  and  HNO2  formed  in  the  heated  channel  rained 
out.  For  extremely  soluble  species  such  as  nitric  and  nitrous 
acids,  this  is  a  reasonable  assumption. 

The  remaining  NO  is  chemically  converted  to  NO2,  primarily 
by  the  reactions 

NO  +  O3  —  * NO2  +  O2 
NO  +  RO2 - *-N02  +  RO 

where  RO2  is  HO2 ,  CH3O2  and  higher  organic  peroxy  radials.  The 
chemical  conversion  occurs  in  2  to  6  minutes  for  temperatures  be¬ 
tween  300  and  240  K  or  altitudes  between  0  and  10  km.  At  somewhat 
higher  temperatures,  before  the  mixing  and  dilution  is  complete, 
the  conversion  occurs  at  a  faster  rate.  During  the  daytime, 
photochemical  steady  state  exists  for  NO  and  NO2 .  The  NO  lost 
chemically  by  reactions  4.9  and  4.10  is  rapidly  regenerated  by 
NO2  photolysis  which  in  turn  regenerates  ozone 

N02  +  hv - ►NO  +  0  \<  420  nm  (4.11) 

O  +  02  +  M - ^03  +  M  (4.12) 

Logan  et  al.  (1981)  show  the  diurnal  variation  of  NOx  and 
HOx  calculated  with  a  1-D  model  for  the  tropics  (15°  N)  and 
midlatitudes  (45°  N)  at  heights  of  0  and  6  km.  The  ground 
level  midlatitude  results  reproduced  in  Figure  8  show  that 
both  the  NOx  and  HOx  chemical  systems  respond  immediately  to  the 
sudden  solar  intensity  changes  at  sunrise  and  sunset.  Clean 
air  measurements  at  Niwot  Ridge,  CO,  have  shown  NO2/NO  ratios 
of  100  ppt/20  ppt  in  the  morning  and  300  ppt/50  ppt  in  the 


(4.9) 

(4.10) 
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FIGURE  8.  Diurnal  behavior  of  the  N0X  AND  HOx  chemical 
systems  in  the  troposphere  at  15°  N  latitude 
equinox,  calculated  with  1-D  model.  The 
upper  set  of  figures  is  for  ground  level  con¬ 
ditions;  the  lower  set  is  for  an  altitude  of 
6  km. 

Source:  Logan  et  al.  1981. 
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late  afternoon  (Kelly  et  al.  1980).  The  rate  of  oxidation  of 
NO  to  NO2#  and  thus  the  ratio  of  NO2/NO,  depends  on  the  local 
NO,  O3  and  organic  peroxy  radical  concentrations,  the  tempera¬ 
ture,  the  sunlight  intensity,  as  well  as  any  local  turbulence 
which  disperses  the  NOx.  A  large  injection  of  NO  can  deplete 
the  local  O3  concentration,  slowing  the  conversion  time  of  NO 
to  NO2  to  several  hours  (Jordan  and  Broderick  1979).  About 
70  percent  of  the  solar  radiation  in  the  visible  wavelength 
region,  where  NO2  absorbs  strongly,  is  reflected  by  dense  cloud 
tops.  During  a  thunderstorm,  turbulence  increases  the  oxidation 
rate  to  NO2  and  the  dark  overhead  sky  reduces  the  NO2  photolysis 
rate,  so-  that  the  NO2  concentration  can  build  up. 

Noxon's  (1976)  measurement  of  NO2  produced  by  lightning 
can  be  considered  a  lower-limit,  estimate  of  the  actual  NOx  pro¬ 
duced,  since  total  conversion  to  NO2  was  assumed  and  no  conver¬ 
sion  to  HNO3  was  considered.  The  N0X  that  is  immediately 
converted  to  HNO3  after  the  lightning  flash  and  is  rained  out, 
is  lost  to  the  atmospheric  NOx  inventory.  It  is  not  necessary 
to  estimate  this  quantity  since  we  are  concerned  with  gas  phase 
NOx  produced  by  lightning  that  can  remain  after  the  thunderstorm 
to  add  to  the  tropospheric  NOx  inventory.  Any  losses  of  NOx 
that  occur  during  the  thunderstorm  but  after  Noxon’s  measure¬ 
ments  do  represent  losses  that  need  to  be  estimated. 

2.  NOx  Loss:  HNO^  Formation  by  NO?  +  OH  +  M 

To  determine  the  efficiency  of  the  gas  phase  reaction  that 
produces  HNO3  it  is  necessary  to  estimate  the  OH  concentration. 
Experiments  by  Campbell  et  al.  (1979),  in  which  chemical  oxida¬ 
tion  of  radioactive ly  labeled  CO  were  used,  indicate  that  under 
daytime  conditions  the  OH  concentration  can  vary  between  3  x  105 
to  3  x  10®  cm“3  in  the  unpolluted  tropospheric  boundary  layer. 
Tropospheric  measurements  by  Wang  and  Davis  (1974),  Wang  et  al. 
(1975),  and  Davis  et  al.  (1976,  1979)  using  laser  fluorescence 
indicate  OH  concentrations  that  are  an  order  of  magnitude 
higher,  but  it  is  possible  that  an  experimental  artifact  caused 
unnaturally  high  values  to  be  measured  (Ortgies  et  al.  1980). 
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Most  tropospheric  OH  is  initially  produced  by  the  mechanism 

03  +  h,  (X  <  320  nm) - ►0(1D)  +  02  (4.13) 

0(1D)  +  H20  - - ►  2  OH  (4.14) 

so  that  one  can  expect  a  variation  in  solar  intensity  to  produce 
a  variation  in  OH  concentration  (Levy  II,  1972).  This  is  corrob¬ 
orated  by  the  ozone  photodissociation  constant  measurements  by 
Dickerson  et  al.  (1979)  (see  Figure  9).  The  variation  of  Of^-D) 
production  with  solar  zenith  angle  was  what  one  would  expect 
from  atmospheric  path  length  considerations.  Superimposed  on 
the  data  were  rapid  spikes  which  corresponded  to,  at  most,  a  50 
percent  reduction  in  O(^D)  production  caused  by  passing  clouds. 
The  UV  intensity  that  produces  the  residual  CK^-D)  is  most 
likely  due  to  scattering  from  the  blue  sky  surrounding  the 
cloud.  An  experiment  in  which  the  solar  disk  was  blocked  out 
but  the  remaining  daylight  was  allowed  to  enter  the  reaction 
vessel  supported  this  hypothesis  (Dickerson  1980). 


TRUE  SOLAR  TIME 

FIGURE  9.  Curve  A,  Data  for  O3  +  h„ - ►  O(^’D)  +  02  photo¬ 

dissociation  constant,  showing  the  effect  of 
passing  clouds.  o  was  calculated  using  atmos¬ 
pheric  absorption  and  Rayleigh  scattering. 
Curve  B  is  exp  (-1/cos  0  ). 

Source:  Dickerson  et  al.  1979. 
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In  order  to  make  estimates  for  the  OH  concentration  that 
exists  during  a  thunderstorm,  one  needs  to  consider  the  amount 
of  light  that  produces  0(1d),  i.e.,  the  295-320  nm  wavelength 
region,  that  can  reach  a  parcel  of  air  under  the  tnundercloud . 
The  optical  density  for  UV  and  near-UV  radiation  is  sufficiently 
high  within  the  cloud  so  that  the  cloud  can  be  considered  opaque 
Scattering  from  the  surrounding  sky  to  the  parcel  of  air  under¬ 
neath  the  cloud  can  be  considered  negligible,  since  most  thunder 
clouds  hang  low  near  the  ground  and  extend  far  vertically,  i.e., 
the  base  of  a  cumulonimbus  cloud  is  typically  0-3  km  above  the 
ground,  and  the  vertical  extent  is  5  km.  Direct  illumination 
beneath  the  cloud  for  the  maximum  3  km  height  cloud  is  only 
possible  for  zenith  angles  greater  than  40°,  which  occurs  at 
local  solar  -times  between  3  p.m.  and  sunset.  Thunderstorms 
generally  occur  in  the  late  afternoon  when  the  Of^D)  production 
rate  is  rapidly  decreasing,  i.e.,  during  the  summer  at  3  p.m. 
solar  time  the  0(1d)  production  rate  was  lower  by  60  percent 
when  compared  to  the  noontime  value  (see  Figure  9).  The  loss 
of  the  photolytic  source  for  OH  in  most  of  the  troposphere 
during  a  thunderstorm  results  in  an  extremely  low  OH  concentra¬ 
tion  if  the  chemical  sources  are  snail. 

During  the  daytime  the  rate  of  interconversion  of  OH  and 
HO2  is  exceedingly  fast  so  that  the  chemical  lifetime  of  OH  is 
dominated  by  the  reactions  that  determine  the  chemical  lifetime 
of  HOx.  The  photolysis  reaction  sequence  4.13  -  4.14  is  the 
main  source  of  tropospheric  HOx .  The  sinks  of  HOx  are  radical 
termination  reactions  such  as 

OH  +  H02 - -H20  +  02  (4.15) 

OH  +  N02  +  M— *-HN03  +  M  aerosol.'*  Products  (4.16) 

rain 


H02  +  H02 


H202  +  02 


aerosol 


products 


(4.17) 


There  are  many  chemical  reactions  or  reaction  sequences  that 
cycle  H0X;  for  example,  the  reactions  of  OH  with  O3,  CO,  CH4, 
and  H2O2  an<^  fc^e  reactions  of  HO2  with  O3  and  NO.  These  proc¬ 
esses  establish  a  steady-state  concentration  of  OH  on  the 
order  of  10^  to  10^  molecules/cm^ . 

The  daylight  ambient  concentrations  of  HO2  and  NO2  are  gen¬ 
erally  about  100  times  greater  than  OH  and  assure  that  the  OH 
is  rapidly  removed  after  the  photolytic  source  of  OH  ceases. 

In  the  absence  of  sunlight  the  HOx  source  term  vanishes  but 
the  chemical  removal  reactions  remain  operative.  This  implies 
an  immediate  drop  in  OH  concentration  and  a  slower  decay  of 
HO2  which  is  verified  by  model  calculations  through  a  solar 
eclipse  (Herman  1979a)  and  diurnal  model  calculations  at  sunset 
(Logan  et  al.  1981),  an  example  of  which  is  shown  in  Fig.  8. 

In  a  darkened  sky  the  OH  concentration  is  not  high  enough  for 
the  NO2  +  OH  +  M  reaction  to  proceed  at  an  appreciable  rate 
during  the  lifetime  of  the  thunderstorm. 

3 .  NOx  Loss.-  HNO 3  Formation  via  N 2O 5 

In  the  absence  of  sunlight,  most  of  the  atmospheric  NO  is 
converted  to  NO2  and  NO3  and  an  equilibrium  between  NO2 ,  NO3  and 
N2O5  is  established.  Experiments  by  Noxon  et  al.  (1978)  show 
that  in  the  stratosphere  this  process  occurs  in  about  40  minutes. 
To  some  extent,  one  would  expect  this  to  also  occur  in  the  tro¬ 
posphere  during  a  thunderstorm  when  the  overhead  sky  is  dark. 

The  N2O5  can  react  with  water  vapor,  present  at  high  concentra¬ 
tions  during  a  thunderstorm,  to  form  HNO3 .  (See  reactions  4.3 
to  4.5.)  The  rate  constant  for  N2O5  +  H2O  is  known  only  as  an 
upper  limit  and  this  discussion  presents  the  maximum  possible 
NO2  loss  to  HNO3  via  N2O5 .  A  further  complication  exists  since 
it  is  possible  that  this  reaction  is  actually  heterogeneous, 
therefore,  aerosols  and  particulates  may  be  necessary  for  HNO3 
formation . 
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The  lifetime  of  NO2  can  be  estimated  from  the  ratio  of  NO2 
to  its  loss  rate  by  this  mechanism 

[N02]  (4.18) 

k5  [N205]  [h2o] 

The  N2O5  concentration  has  not  been  measured  in  the  troposphere 
but  it  can  be  inferred  from  NO2  and  NO3  concentrations  and  the 
equilibrium  constant  which  have  been  measured.  An  upper  limit 
for  the  tropospheric  vertical  column  density  of  N03  of  0.2  x  1014 
cm“2  has  been  measured  by  Noxon  et  al.  (1978).  In  this  experiment 
NO2  concentrations  less  than  1  ppb  were  required  to  be  consistent 
with  calculations.  In  a  similarly  clean  atmosphere  Noxon  (1978) 
measured  N02  column  densities  of  1015  cm"2  in  the  troposphere, 
which  corresponds  to  0.2  ppbv  at  the  surface.  The  equilibrium 
constant  for  reaction  4.4  shows  an  extreme  temperature  dependence, 
i.e.,  the  calculated  N2O5  concentration  increases  by  a  factor 
of  10  for  a  20-degree  temperature  decrease  of  320  to  300  K 
(Connell  and  Johnston  1979).  For  this  same  temperature  range 
the  water  vapor  pressure  drops  by  a  factor  of  3.  With  a  typical 
midlatitude  water  vapor  content  during  a  storm  of  about  3  x  1016 
molecules/cm3  the  lifetime  from  equation  4.18  is  about  500  hrs  at 
300  K.  Using  the  most  extreme  water  vapor  content  recorded  at 
ground  level  for  the  tropical  regions  of  3  x  1037  molecules/cm3 
reduces  all  the  lifetime  estimates  by  a  factor  of  10.  At  minimum 
tropospheric  temperatures,  i.e.,  between  10  to  20  km,  N2O5  produc¬ 
tion  becomes  more  favorable;  however,  the  H2O  vapor  concentration 
is  lower  by  four  orders  of  magnitude.  The  resulting  removal  of 
NO2  via  N2O5  is  minimal  high  in  the  troposphere  if  the  N2O5 
concentration  has  not  increased  more  than  3  orders  of  magnitude. 
Such  as  increase  is  unlikely,  since  the  local  NO2  concentration 
sets  an  upper  limit  on  the  N2O5  concentration  present  at  low 
temperatures.  Also,  the  arguments  which  depend  on  an  equilibrium 
between  N2O5,  NC^*  and  NO3  are  uncertain  since  one  cannot  assume 
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an  equilibrium  exists  in  the  cold  troposphere  where  the  rate  of 
thermal  decomposition  of  N2O5  is  very  slow. 

The  best  estimate  that  can  be  made  for  the  conversion  of 
NO2  to  HNO3  in  the  absence  of  sunlight  is  that  this  conversion 
is  minimal  throughout  the  troposphere  during  the  lifetime  of  a 
thunderstorm,  but  it  is  difficult  to  quantify  it  more  exactly. 
It  should  be  possible  to  make  a  better  assessment  when  the 
tropospheric  N2O3  concentration  has  been  measured  directly. 
Conclusions  based  on  NO3  and  NO2  data  alone  are  difficult  to 
use  due  to  the  sensitivity  of  the  NO,  NO2 ,  NO3 ,  and  N2O5  chemi¬ 
cal  system  to  the  vertical  temperature  profile. 

4 .  N0X  Loss:  Dry  Deposition 

The  fourth  mechanism  for  the  loss  of  NO  and  NO2  created 
by  lightning  is  dry  deposition  onto  ground  or  water  surfaces. 
The  deposition  velocity  which  is  used  to  quantify  the  rate  of 
mass  transfer  from  the  gas  phase  to  the  surface  is  typically 
between  0.1  and  1  cm/s  at  a  reference  height  of  1  km  or 
at  the  top  of  the  logarithmic  wind  profile  within  the  planetary 
boundary  layer.  The  value  of  the  deposition  velocity  depends 
on  a  number  of  factors;  for  example,  the  altitude,  the  exact 
nature  of  the  interaction  of  the  molecule  with  the  surface, 
the  temperatjre,  the  wind  speed,  etc.  The  following  deposi¬ 
tion  velocities  for  species  of  interest  in  the  NOx  chemical 
system  have  been  recently  measured  over  a  grassland:  HNO3  = 

0.34  +  0.25,  03  =  0.67  +  0.15,  NOx  =  0.73  +  0.24  cm/s  at  a 
reference  height  of  1  m  (Kasting,  1980).  For  a  1-km  mixing 
heiqht  of  the  daytime  boundary  layer,  a  range  for  the  deposi¬ 
tion  velocity  of  0.1  to  l  cm/s  results  in  a  deposition  rate 
of  10-5  to  10"^  seconds"!  or  a  i/e  lifetime  between  1-10 
days.  For  a  chemical  species  as  reactive  as  NO2 ,  the  faster 
deposition  rate  is  probably  more  appropriate. 

5.  NOx  Loss:  Aerosol  Scavenging  and  Rainout 

The  final  mechanism  for  loss  of  NO2  is  heteroqeneous  re¬ 
actions  on  particles  or  water  droplets.  It  is  unlikely  that 


significant  amounts  of  gaseous  NO2  will  dissolve  into  water 
droplets  to  be  rained  out  during  the  lifetime  of  the  thunder¬ 
storm,  since  the  various  factors  governing  the  solubility  of 
gaseous  NO2  are  unfavorable.  Dissolution  involves  three  equi¬ 
librium  processes:  (1)  partitioning  between  the  gaseous  and 
aqueous  phases  which  is  parameterized  by  the  Henry's  Law  co¬ 
efficient,  (2)  an  aqueous  phase  equilibrium  of  the  NO2  to  form 
chemically  associated  oxyacids  such  as  HNO2  and  HNO3,  and  (3) 
an  acid  dissociation  equilibrium  in  which  the  ionic  species 
H+,  NC>2~  or  NO3-  are  formed.  The  equilibrium  concentration 
calculated  by  Henry's  Law  indicates  that  the  gaseous  NO2  can 
enter  the  liquid  phase  to  form  a  0.007  molar  solution  per  at¬ 
mosphere  NO2  (Lee  and  Schwartz  1981a).  Partitioning  into  the 
atmosphere  is  favored;  the  NO2  experiences  liquid  phase  resis¬ 
tance  to  mass  transfer  (Eisenreich  et  al.  1981).  The  extent  to 
which  the  NO2  forms  dissociated  acid  within  the  liquid  depends 
on  the  kinetics  of  the  following  reaction. 

2  N02(aq)  +  H20(aq) - «-2H+  +  N02"  +  N03~  (4.19) 

Lee  and  Schwartz's  (1981a)  recent  kinetic  study,  done  with  NO2 
partial  pressures  comparable  to  those  found  in  the  polluted 
troposphere,  showed  that  the  time  constant  for  the  above  re¬ 
action  is  at  least  ten  times  slower  than  the  time  for  convec¬ 
tive  or  diffusive  mass  transport  within  the  liquid.  This  in¬ 
dicates  that  NO2  dissolution  is  limited  by  both  the  rate  of 
physical  transfer  from  the  gas  phase  to  the  aqueous  phase  and 
by  the  rate  of  ion  formation  in  solution  to  a  sufficient  extent 
so  that  very  little  of  the  NO2  produced  during  the  lightning 
storm  is  immediately  dissolved  and  rained  out. 

Recently  an  experimental  study  was  presented  in  which  the 
chemistry  of  NO2  with  realistically  simulated  clouds  was  studied. 
(Miller  and  Gertler  1981).  For  0.1  ppm  NO2  in  clean  air,  the 
rate  of  conversion  to  NO^-  in  the  cloud  was  0.04  percent/hr. 

This  corresponds  to  a  lifetime  of  1700  hrs  for  NO2. 
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It  is  also  possible  that  the  dissolution  of  NCV)  may  be 
aided  by  the  presence  of  oxidizing  or  reducing  agents  present 
in  the  rain  droplets  or  on  the  surface  of  particles.  Although 
heterogeneous  chemistry  :1s  not  well  understood,  some  progress 
has  been  made  by  the  environmental  requirement  to  remove  N0X 
from  power  plant  and  industrial  plant  effluents.  Most  N0X 
removal  schemes  include  the  use  of  a  catalyst  or  a  chemical 
reaction  sequence  which  first  oxidizes  or  reduces  the  N0X  to 
enhance  the  solubility  into  the  liquid  phase  where  it  is  even¬ 
tually  converted  to  a  gas  or  to  an  insoluble  solid  which  pre¬ 
cipitates  (Faucett  1977).  Studies  by  Lee  and  Schwartz  (1981b) 
indicate  that  the  uptake  of  NO2  by  liquid  water  is  potentially 
important  for  catalyst  concentrations  on  the  order  of  10 “7  M 
and  a  fast  rate  constant  for  reaction  4.19.  They  found  that  the 
Fe  catalyzed  reaction  is  not  important  but  the  reaction  of  NO2 
with  dissolved  S(IV)  is  potentially  important.  It  is  difficult 
to  estimate  to  what  extent  such  heterogeneous  chemistry  enhances 
the  dissolution  of  NO2  into  raindrops  in  the  troposphere,  since 
the  concentration  of  catalysts  and  othe'  dissolved  species  in 
raindrops  is  uncertain. 

B.  HNO3  LOSS:  RAINOUT ,  DRY  DEPOSITION  AND  HETEROGENEOUS 
SCAVENGING 

The  small  amount  of  the  gaseous  HNO3  that  is  formed  during 
the  lightning  storm,  as  well  as  the  residual  HNO3  formed  before 
the  storm,  is  quickly  rained  out.  Table  12  summarizes  the 
partitioning  between  gaseous  and  aqueous  HNO3  calculated  from 
equilibrium  considerations.  The  aqueous  phase  is  stronqly 
favored.  Table  13  shows  that  the  scavenging  coefficient  for 
HNO3  into  cumulus  clouds  is  fast  enough  that  within  a  cloud, 
gaseous  HNO3  is  quickly  dissolved  within  water  droplets  (Levine 
and  Schwartz  1981).  The  scavenging  coefficients  for  a  heavy 
rainfall  are  also  fast  enough  that  any  residual  HNO3  vapor  is 
quickly  picked  up  by  the  falling  water  droplets  and  rained  out. 


TABLE  12.  EQUILIBRIUM  OF  HNO3  WITH  WATER* 


pHN03  (atm) 

(HN03]g**  mol  cm-3 

[HN03]aq/[HN03]g 

10“8 

2.6  x  1011 

50 

10-10 

2.6  x  109 

500 

Source:  Levine 

and  Schwartz  1981. 

*The  water  vapor  content  assumed  is  10-6  g  cm-3. 

**Note  that  the  measured  fHN03]g  *-n  the  troposphere  range 
from  0.1  -  5  ppb  (2.6  x  1G9  -  1.3  x  1011  mol  cm-3) 
(Kelley  et  al.  1979). 


TABLE  13.  HNO3  SCAVENGING  BY  CLOUDS  AND  RAINFALL 


(scavenging  Medium 

Scavenging  Coefficient 

Scavenging  Time 

Cumu'us  clouds 

0.2  s~l 

5  seconds 

Light  rain 

1  mm  hr~l 

1.3  x  10-5  s*1 

21  hours 

Heavy  rain 

25  mm  hr-3 

1.5  x  IO-3  s-1 

11  minutes 

Long-term  average: 
includes  rain  and 
dry  deposition. 

(1-8)  x  10"6  s'1 

1-11  days 

Aerosol  Scavenging  by  Rainfall 

Medium  rain 

10  -  18  mm  hr-l 

1  -  3.3  x  10~3  s"1 

5-15  minutes 

Source:  HNO3  data  compiled  from  Levine  and  Schwartz  1981. 

Aerosol  data  from  Radke  et  al.  1980. 

■  -  — ■  ■  -  -  ~  ■"  ’  mmm‘ 

The  importance  of  heterogeneous  rainout  was  shown  by  a  tropo¬ 
spheric  modelling  calculation  for  HNO3  by  Fishman  and  Crutzen 
(1977).  An  order  of  magnitude  decrease  in  tropospheric  HNO3 
concentration  was  noted  when  a  heterogeneous  removal  process 
with  a  rate  of  2  x  10“®  s“l  was  added  for  the  lowest  6  km.  The 
calculation  by  Levine  and  Schwartz  (1981)  shows  that  the  long¬ 
term  average  for  scavenging  by  falling  raindrops  and  dry  depo¬ 
sition  is  (1  to  8)  x  10"6  s“l. 

In  a  heavy  rainstorm  the  scavenging  efficiency  of  gaseous 
HNO3  by  falling  raindrops  and  the  scavenging  efficiency  of 
aerosol  particles  by  rain  are  comparable.  This  is  shown  in 
Table  13.  In  addition,  if  one  assumes  that  deposition  of  HNO3 
onto  atmospheric  aerosol  particles  occurs  at  the  same  rate  as 
dry  deposition  to  the  ground,  the  removal  of  HNO3  vapor  by 
aerosols  is  very  effective. 

Overall,  the  removal  of  HNO3  by  various  heterogeneous  proc¬ 
esses  is  very  efficient.  This  is  what  one  would  expect  since 
20  to  40  percent  of  the  acid  content  in  rain  over  polluted 
areas  is  HNO3  (Likens  et  al.  1979).  Also,  measurements  of 
the  HNO3  content  in  air  before  and  after  rainfall  shows  that 
rainfall  is  a  very  effective  removal  process.  Previous  workers 
have  tried  to  correlate  the  HNO3  measured  in  rain  with  lightning 
activity  with  no  success  (Viemeister  1960;  Visser  1960;  Reiter 
1970).  This  is  what  one  would  expect  if  the  NOx  produced  by 
lightning  is  not  converted  to  HNO3  during  the  lifetime  of  the 
thunderstorm. 
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V.  UNCERTAINTIES  AND  RECOMMENDATIONS 


A.  UNCERTAINTIES 

There  are  several  major  sources  of  uncertainty  apparent  r.-n 
the  estimate  of  the  lightning  contribution  to  the  global  N0X 
budget: 

•  The  NOx  production  efficiency  for  both  cloud- to-ground 
lightning  and  intracloud  discharges. 

•  The  injection  height  of  N0X  produced  during  a  light¬ 
ning  storm,  and  its  vertical  distribution. 

•  The  normalization  of  the  global  lightning  rate  meas¬ 
ured  by  satellite,  especially  in  the  tropics. 

•  The  ratio  of  cloud  discharges  to  ground  flashes,  also 
with  particular  emphasis  in  the  tropics. 

Noxon  (1976,  1978)  estimated  the  uncertainty  of  the  NO2  produc¬ 
tion  rate  to  be  within  one  order  of  magnitude.  This  does  not 
include  the  systematic  error,  if  any,  introduced  by  extrapolating 
Noxon's  midlatitude  measurement  to  the  tropics,  and  by  assuming 
that  the  HO2  measured  by  Noxon  is  equivalent  to  N0X. 

It  should  be  emphasized  here  that  in  some  cases  assumptions 
have  been  made  because  direct  experimental  evidence  does  not 
exist.  For  example,  the  estimate  of  N0X  production  by  intra¬ 
cloud  discharges  has  been  based  on  the  argument  that  the  amount 
of  NOx  produced  is  directly  proportional  to  the  discharge  energy 
and  that,  on  the  average,  cloud  discharges  are  about  one  tenth 
as  energetic  as  the  ground  flash.  The  mean  stroke  length  for 
intracloud  and  cloud- to-ground  lightning  was  assumed  to  be  equiva¬ 
lent.  Even  though  some  exceptionally  long  intercloud  discharges 
have  been  observed,  e.g. ,  50  km,  much  cloud  lightning  is  of  the 
the  shorter,  intracloud  type.  The  maximum  height  of  NOx  injection 
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was  assumed  to  follow  the  height  of  local  tropopause,  but  again, 
no  experimental  evidence  exists  to  support  this  conclusion. 

\  A  large  uncertainty  is  introduced  by  current  normalization 

techniques  employed  for  satellite  experiments.  Simultaneous 
satellite  and  ground-based  measurements  provide  good  calibra¬ 
tion  points,  provided  that  some  tropical  as  well  as  midlatitude 
I  coverage  is  included.  For  accurate  global  counting,  the  effi¬ 

ciency  of  intracloud  detection  needs  to  be  carefully  considered, 
since  certain  satellite  and  ground-based  techniques  selectively 
discriminate  against  the  less  energetic  events  such  as  intracloud 
discharges. 

Another  uncertainty  is  introduced  through  the  large  scatter 
of  the  Nc/Ng  ratio  evident  in  Figure  1.  In  part,  this  is  due 
to  an  inadequate  statistical  sample  and  in  part  to  natural  varia¬ 
tions  over  different  geographical  terrain.  The  deficiency  can 
be  remedied  with  more  satellite  experiments  properly  equipped 
with  sensors  equally  sensitive  to  cloud  discharges  and  ground 
flashes.  Special  attention  should  be  given  to  the  tropical  re¬ 
gions  where  current  data  is  sparse  but  where  most  of  the  light¬ 
ning  activity  occurs. 

The  overall  uncertainty  has  been  estimated  to  be  approxi¬ 
mately  one  order  of  magnitude,  reflecting  the  accumulation  of 
uncertainty  in  the  measured  quantities  only.  For  quantities 
where  no  measurements  exist,  such  as  the  N0X  production  effi¬ 
ciency  by  cloud  discharges  and  the  injection  height,  it  is 
not  possible  to  estimate  the  uncertainty. 

B.  RECOMMENDATIONS 

Most  of  the  uncertainties  mentioned  can  be  resolved  using 
previously  developed  experimental  systems.  For  example,  a  set 
of  experiments  capable  of  measuring  in-situ  N0X  at  various  al¬ 
titudes  during  a  thunderstorm  could  eliminate  the  uncertainty 
of  the  injection  height  and  the  uncertainty  of  the  N0X  produc¬ 
tion  efficiency  by  cloud  discharges.  Similar  experiments  are 
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routinely  made  for  environmental  monitoring  of  N0X  in  ambient 
air  at  ground  level.  The  NO  concentration  has  been  measured 
from  an  aircraft  during  the  GAMETAG  experiments  (Davis  1980) 
and  NO  and  NO2  have  been  measured  several  times  in  balloons 
(see  for  example  McFarland  et  al.  1979). 

The  following  recommendations  are  made  here: 

•  The  F-106B  aircraft  used  in  the  NASA  Langley  Storm 
Hazard  Project,  which  is  currently  equipped  to  fly 
within  thunderstorms,  could  be  easily  adapted  to 
accommodate  an  N0X  monitoring  system.  Several 
measurements  should  be  made  at  various  altitudes 
within  the  thundercloud,  at  its  perimeter,  and 
outside  the  cloud  in  order  to  sample  regions  of  net 
updrafts  and  downdrafts.  This  single  experiment 
would  answer  both  the  question  of  NOx  production 
efficiency  by  cloud  discharges  as  well  as  the  net 
injection  height  during  and  after  the  thunderstorm. 

A  global  aircraft  study  is  not  needed  at  this  time; 
global  behavior  could  be  inferred  from  other 
parameters. 

•  The  NOx  measurements,  such  as  Noxon's  column  densities, 
are  needed  in  the  tropics  to  determine  whether  the 
production  efficiency  by  lightning  is  different  for 
tropical  thunderstorms. 

•  Experimental  observation  of  lightning  by  satellite 
presents  a  comprehensive  and  consistent  global  dis¬ 
tribution,  but  a  consensus  for  the  globally  averaged 
lightning  rate  does  not  exist.  The  problem  of  cali¬ 
bration  of  the  satellite  lightning  measurements  should 
be  specifically  addressed.  Satellite  experiments 
should  be  combined  with  ground-based  measurements  to 
provide  an  absolute  calibration  for  the  flash  rate 

in  both  the  tropics  and  the  midlatitudes.  Preferably, 
experimental  systems  capable  of  detecting  both  cloud 
discharges  and  ground  flashes  should  be  chosen. 


•  Several  measurements  oE  the  ratio  of  cloud  discharges 
to  ground  flashes  are  needed  in  the  tropics  to  reduce 
the  scatter  evident  in  currently  available  data.  In 
general,  tropical  regions  are  poorly  represented  by 
individual  experimental  measurements,  but  it  should 
be  recognized  that  most  lightning  occurs  there. 

An  experiment  in  which  all  NOy  species,  i.e.,  NO,  NO2 , 
HNO3,  HN02»  and  possibly  NH3  and  NO3,  will  be  measured  simul¬ 
taneously  during  a  thunderstorm  is  currently  being  prepared 
at  the  Oeorgia  Institute  of  Technology  (M.  0.  Rodgers,  private 
communication).  The  concentrations  will  be  monitored  in  a 
parcel  of  air  before,  during,  and  after  a  thunderstorm,  using 
laser-induced  fluorescence  combined  with  photofragmentation 
techniques  when  necessary  (Rodgers  et  al.  1980).  Although 
this  experimental  system  cannot  be  routinely  employed  or  placed 
in  an  aircraft  because  it  is  too  expensive,  bulky,  and  heavy, 
several  individual  experiments  would  answer  questions  about 
chemical  production  efficiencies  by  lightning,  the  partitioning 
among  the  NOy  species  in  the  thunderstorm  environment,  and  the 
effectiveness  of  rainout  for  the  various  chemical  species.  The 
direct  concentration  measurements  to  be  obtained  from  this 
experiment  complement  Noxon's  measurements  of  the  NO2  column 
density. 
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